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Abstract 

In these notes we review the material presented at the summer school on “Mathematical 
Physics, Analysis and Stochastics” held at the University of Heidelberg in July 2014. We con¬ 
sider the time-evolution of quantum systems and in particular the rigorous derivation of effective 
equations approximating the many-body Schrddinger dynamics in certain physically interesting 
regimes. 


1 Introduction 

Systems of interest in physics and other natural sciences can be described at the microscopic and the 
macroscopic level. Microscopically, a system is described in terms of its elementary constituents and 
their fundamental interactions. While such a description is very accurate, it is typically not well- 
suited for computations because of the large number of degrees of freedom. Examples of microscopic 
theories include Newton’s theory of classical mechanics, Schrodinger’s quantum mechanics, quantum 
electrodynamics and Einstein’s general relativit50. On the other hand, a macroscopic description 
of the system does not resolve the constituents and only takes into account effective interactions. 
It focuses on macroscopically observable quantities which arise from the collective behavior of the 
system and are of interest for the observer. Such a description is less accurate but it is much more 
accessible to computations. Examples of macroscopic theories are Boltzmann’s kinetic theory of 
gases, the Navier-Stokes and the Euler equations of hydrodynamics, the Hartree and Hartree-Fock 
theory, the BCS theory of superconductors and superfluids, the Ginzburg-Landau theory, the Gross- 
Pitaevskii theory of Bose-Einstein condensation and the Vlasov theory of plasma physics. 

Because of the great importance of effective macroscopic theories for making qualitative and 
quantitative predictions about the behavior of physically interesting systems, a key goal of statistical 
mechanics is to understand their emergence from microscopic theories in appropriate scaling regimes 

^Of course we are not claiming these theories to be absolutely fundamental from the view of a physicist. It would 
be more correct to consider them as different levels between fundamental and effective, and which theory we call 
effective and which fundamental depends on the pair we are looking at. For example, we could also consider Newtonian 
mechanics as a macroscopic theory arising as an effective theory from quantum mechanics. On the next level we could 
view nonrelativistic quantum mechanics as an effective theory arising from the Standard Model. 
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(also called limits, even though we often think of the parameter as being large but finite). Here 
mathematical physics can and should play a central role to put the effective theories, which are often 
obtained merely by heuristic and phenomenological arguments, on solid grounds and understand the 
range and the limits of their validity. 

Let us present two examples of physical systems which can be described by effective equations 
that can be rigorously derived from microscopic theories in appropriate limits. 

Large atoms and molecules. We consider a quantum-mechanical system of N electrons and M 
nuclei of charges Zi,..., Zm > 0 located at positions i2i,..., Rm £ We assume the system to 
be neutral, i. e. N = simplicity we work in the Born-Oppenheimer approximation, i. e. 

we keep the nuclei fixed (only the electrons are dynamical particles in this approximation). At zero 
temperature the system is in its ground state, with energy 


E{N) 


min 




where Hjsf denotes the Hamilton operator 


N 


Hn = Y1 
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ZiZj 

\Ri — Rj 


( 1 . 1 ) 


( 1 . 2 ) 


Notice that acts on the subspace L^(M^'^) of L^(M^'^) consisting of all functions that are anti¬ 
symmetric with respect to permutations of the N electrons. (Of course E{N) and H]\f also depend 
on M, on the charges Zi,, Zm and on the positions Ri,..., Rm)- Observe that the last term 
on the r. h. s. of (ll.2p is just a constant representing the interaction among the nuclei. Already for 
~ 20 it is extremely difficult to compute the ground state energy E{N) numerically since the 
eigenvalue equation one has to solve is a partial differential equation in 3N coupled variables. 

Thomas and Fermi postulated already in the early stages of quantum mechanics at the end of 
the 1920s that the ground state energy E{N) can be approximated by 


E(N) ~ ETpiN) = inf STpip) 

p>0,||p||i=7V 


(1.3) 


with the Thomas-Fermi functional 


^Tp{p) = pCTP 

5 


M 

p^^^{x)dx — Zi 
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p{x)p{y) 
|x - y\ 


M 

dxdy + 

i<j 


ZiZj 

\Ri — Rj 


Notice that on the r. h. s. of ()1.3I1 we are looking for a function p G L^(M^); as a consequence, in 
terms of numerical computations, the minimization problem (II.3p is much simpler than the original 
problem (jl.ip despite the fact that the resulting Euler-Lagrange equations are nonlinear. 

In [49j Lieb and Simon proved that the approximation (II.3p becomes exact in the limit of large 
N. More precisely, they showed that, while E{N) and Etp{N) are both proportional to 

\E{N) - Etp{N)\ < CAf7/3-i/30 


for an appropriate constant C > 0. This gives a mathematically rigorous derivation of the Thomas- 
Fermi theory, and it tells us how big N should be in order for Etp{N) to be a good approximation 
of the true quantum energy (later, better bounds have been obtained [HI [58l [59]: one knows that 
the error with respect to Thomas-Fermi theory is of the order in the limit of large N). 
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Kinetic theory of dilute gases. Consider now a gas of N classical particles moving according to 
Newton’s equations 


Xj{t) = Vj{t) 

N 

vj{t) = - - xj{t)) 


(1.4) 


for j = 1,... ,N. Here H is a short range (compactly supported), regular potential. Eq. (11.41) is a 
system of 6N coupled ordinary differential equations. Given appropriate initial data, it is known to 
have a unique solution for all t G M. However, since the number of particles N is typically extremely 
large, it is almost impossible to deduce from (11.41) interesting qualitative or quantitative properties 
of the solution. 

At the beginning of the twentieth century Boltzmann, based on clever heuristic arguments, pro¬ 
posed to describe the dynamics of the gas by the nonlinear equation 

dtft{x,v) +v - K^ftix^v) = J dv' B{v - v']u;) {ft{x, Vout) ft - ft{x,v)ft{x,v')) (1.5) 

for the phase-space density ft{x, v), which should measure the number of particles at time t that are 
located close to a: G and have velocity close to u G M^. Here 

Vont =V+UJ-{v'- v)uj 
v'out =v'-Uf(v' - v)uj 


are the velocity of two particles emerging from the collision of two particles with velocities v, v' and 
collision vector cj G 5^. Boltzmann’s equation is a partial differential equation in only six variables, 
and is therefore much more accessible to computations than (HH). 

For several years after the work of Boltzmann, his equation was not accepted by the physics 
community. In contrast to Newton’s equations, Boltzmann’s equation m is not time-reversal 
invariant. It took a while to understand that this fact does not contradict the validity of Boltzmann’s 
equation, but only restricts the set of initial data for which ([ES]) is a good approximation. 

In [36] Grad realized that Boltzmann’s equation becomes correct in the low-density limit, where 
the density of particles p is very small, N is very large, and Np'^ is fixed of order one (the low-density 
or Boltzmann-Grad limit). Finally Lanford proved in |45] that indeed in this limit Boltzmann’s 
equation can be rigorously derived from Newton’s equations, at least for sufficiently short times 
(recently |55l|3l] have extended the validity of (|1.5p to a larger class of interaction potentials). 

Plan of the notes. Notice that the first example we discussed above (the Thomas-Fermi theory 
for large atoms and molecules) is based on many-body quantum mechanics while the second one 
(kinetic theory of gases in the low density limit) is based on classical Newtonian mechanics. Observe, 
moreover, that in the first example we were interested in an equilibrium property of the system (its 
ground state energy) while in the second case we considered a non-equilibrium problem (the time- 
evolution of the gas). 

In these notes we will focus on the derivation of time-dependent effective theories (non-equilibrium 
question) approximating many-body quantum dynamics. In the rest of this section, we will briefly 
recall the main properties of many-body quantum systems and their time evolution. In Section 
[2] we will then introduce the mean-field regime for bosonic systems and we will explain how the 
many-body dynamics can be approximated by the Hartree equation in this limit. In Section [3] we 
will present a method, based on the use of coherent states and inspired by [lUl ES], to rigorously 
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prove the convergence towards the Hartree dynamics. The fluctuations around the Hartree equation 
will be considered in Section 01 In Section 0 we will discuss a more subtle regime, in which the 
many-body evolution can be approximated by the nonlinear Gross-Pitaevskii equation. In Section 
El we will discuss fermionic systems (characterized by antisymmetric wave functions). The fermionic 
mean-field regime is naturally linked with a semiclassical regime, and we will prove that the evolution 
of approximate Slater determinants can be approximated by the nonlinear Hartree-Fock equation. 
Finally, in Section [3 we will consider the same fermionic mean-field regime but this time we will 
focus on mixed quasi-free initial data approximating thermal states at positive temperature. 

Wave functions and observables. We will describe quantum systems of N particles in three 
dimensions through a complex-valued wave function ipN ^ ,dxi,... dx^) with ||?/>Ar ||2 = 1- 

The arguments (xi ,... ,xn) £ of the wave function ipN describe the position of the N particles. 
Since it does not play an important role in our analysis, we always neglect the spin of the particles. 

Observables of the quantum system are associated with self-adjoint operators on the Hilbert space 
dxi... dxN)- If A is such an operator, the spectral theorem allows us to write it as 

A = j XdEAiX), 

where dEA{X) is the spectral measure associated with A. Assuming for simplicity that A has purely 
discrete spectrum, we find A = Xj\ifj){ifj\, where Xj are the eigenvalues and ipj the eigenvectors 
of A. The expectation of the observable A in the state ?/>Ar is then given by the inner product 

j 

This expression for the expectation of A leads to the following interpretation: the eigenvalues Xj 
are the possible outcomes of a measurement of A while is the probability that a mea¬ 

surement of A produces the outcome Xj. The positions of the N particles are associated with mul¬ 
tiplication operators; hence \^|^]\^{xl,... ,X]y)\‘^ is the probability density for finding particles close 
to (xi,...,xjv) G . The momenta of the particles are associated with differential operators 
Pj = —iS/xj- Hence, if denotes the Fourier transform oi if n-, |^iv(pij • • • ,Piv)P is the probability 
density for finding particles with momenta close to (pi ,... . The fact that ipN determines 

the probability distribution of all observables of the system is an important feature of quantum 
mechanics (it leads for example to Heisenberg’s uncertainty principle). 

Statistics. In the following we will study systems of N indistinguishable particles. In this case, 
there are important restrictions on the behavior of the wave function with respect to permutations. 
There are two different types of particles. Bosons are characterized by wave functions which are 
symmetric with respect to permutation, i. e. tpNixTr{i)i • • • ja^ 7 r(iv)) = f’Nixi,... ,xn) for all permuta¬ 
tions vr G Sisf. Fermions, on the other hand, are characterized by antisymmetric wave functions, i. e. 

... ,x^(jv)) = (XTrip]\f{xi,... ,xn), where Utt is the sign of the permutation vr G SM In these 
notes, we will consider both the time-evolution of bosonic and fermionic systems; as we will see, the 
different behavior with respect to permutations has important consequences on the dynamics. 

Dynamics. The time evolution of quantum systems is governed by the Schrodinger equation 

idt’4’N,t = HN’ipN,t (1-6) 

^Bosonic and fermionic wave functions are the only one-dimensional irreducible representations of the permutation 
group Sn', an irreducible representation of Sn with dimension larger than one would lead to redundancies in the 
description of the system and is therefore meaningless, from the point of view of physics. 
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for the A^-particle wave function (the subscript t indicates the time dependence of On 

the r. h. s. of (II.6p . H]\f is a self-adjoint operator on ,dxi ...dxjy) known as the Hamilton 

operator (or Hamiltonian) of the system. In these notes, we will consider Hamilton operators with 
two-body interactions, having the form 

N N 

= + ^ext(Xj)) + A ^ V{Xi - Xj) (1.7) 

j=l i<j 

for appropriate functions I4xt (the external potential) and V (the interaction potential) and for a 
coupling constant A G M (which we introduce here for later convenience). The snm of the Laplace 
operators is the kinetic part of the Hamiltonian and generates the evolution of free particles. 

The Schrodinger equation (II.6p is linear and can always be solved by where 

denotes the unitary one-parameter group generated by the self-adjoint operator Hn. This 
implies that the well-posedness of the Schrodinger equation is not an issue (although sometimes it 
can be difficult to prove the self-adjointness of the Hamilton operator). Nevertheless, for N ^ 1, 
it is essentially impossible to extract useful qualitative and qnantitative information from m that 
goes beyond the existence and uniqueness of solutions. According to the philosophy outlined above, 
we look for simpler effective equations which approximate (fLGD in interesting regimes. 

Acknowledgments. We would like to thank Manfred Salmhofer and Christoph Kopper for organiz¬ 
ing the summer school “Mathematical Physics, Analysis and Stochastics” and for encouraging us to 
write up these notes. The work of M. Porta and of B. Schlein has been snpported by the ERC grant 
MAQD-240518. N. Benedikter has been partially supported by the ERC grant CoMBoS-239694 and 
by the ERC Advanced grant 321029. 

2 Mean-field regime for bosonic systems 

One of simplest non-trivial regimes in which it is possible to approximate the many-body dynamics by 
an effective equation is the mean-field limit for bosonic systems. In the mean-field regime, particles 
experience a large number of weak collisions, whose cumulative effect can be approximated by an 
average mean-field potential. To realize the mean-field regime, we consider a system of N bosons 
with a Hamilton operator of the form (II.7p . in the limit of large N and small coupling constant A, 
with NX of order one. This last condition guarantees that the total force on each particle is of order 
one, and therefore comparable with the inertia. In other words, we consider the dynamics generated 
by the mean-field Hamiltonian 

N N 

+ bixt(xj)) + Xj) (2.1) 

j=l i<j 

acting on the Hilbert space the subspace of consisting of permutation symmetric 

functions, in the limit of large N. 

Initial data. The choice of the initial data is dictated by physics. In typical experiments a Bose 
gas is initially trapped by an external confining potential. To study the dynamics of the gas out of 
equilibrium, we consider the reaction of the system to a change of the external field. In other words, 
we are going to consider initial data given by equilibrium states of a Hamiltonian of the form (12.11) . 
with I4xt modeling the external traps. In particular, at zero temperature, we are interested in initial 
data close to the ground state of m- Under appropriate assumptions on the interaction potential 
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V it is known that the ground state of (I2.ip can be approximated, in the limit of large N, by a 
factorized wave function; i. e. for an appropriate (/? E {ip is the minimizer of the 

Hartree energy functional). We are interested, therefore, in the solution of the Schrodinger equation 

idtipN,t = ( 2 . 2 ) 

for approximately factorized initial data V’iv.o — ■ Note that the external potential in the 

operator appearing on the r. h. s. of (| 2 . 2 p is typically different from the external potential in the 
trapping Hamiltonian, whose ground state is approximated by (otherwise, the dynamics would 
be trivial). For example the initial data would be taken as an approximation to the ground state in 
a harmonic trap, and the evolution of this initial data would be studied after switching off the trap, 
i. e. Fext = 0 . 

Hartree equation. Of course, since H]\[ is an interacting Hamiltonian, the dynamics does not 
preserve the factorization of the many-body wave function. Still, since the interaction is weak, we 
can expect factorization to be approximately (in a sense to be specified later) preserved, in the limit 
of large N. In other words, we can expect that for N ^ 1 

N 

1pN,t{xi, . . .,Xn) ^ (2-3) 

i=i 

for an evolved one-particle wave function pt. Assuming ()2.3p . it is easy to derive a self-consistent 
equation for the evolution of the one-particle wave function pt- In fact, factorization of the Wparticle 
wave function means, in probabilistic terms, that the particles are distributed in space according to 
the density independently of each other. The law of large numbers then suggests that the total 
potential experienced, say, by the j-th particle can be approximated by 

^^F(xi-xj)~ f V{xj-y)\pt{y)\‘^dy = {V*\p\‘^){xj). 

Hence, pt must satisfy the Hartree equation 

idtpt = (-A -b Ve^t)pt + (F * \pt?‘)pt (2.4) 

where the many-body interaction has been replaced by the effective one-particle potential F * \pt\‘^, 
making \pt\'^ a nonlinear equation. 

Reduced densities. To explain in which sense we can expect (USD to hold true, we introduce 
the notion of reduced density matrices (also known as reduced densities). The one-particle reduced 
density associated with the wave function 'ipN,t is defined by 

InI = ^Tr2,3,...,Ar|V’Ar,t)(V’Ar,t|, 

where \'tpN,t){'ipN,t\ denotes the orthogonal projection onto 'ipN,t and Tr 2 ,...,Ar is the partial trace over 
the last {N — 1) particles. In other words, the one-particle reduced density 7 ^ ^ is defined as the 
non-negative trace-class operator on L^(M^) with integral kernel 

y) = N j dx2 ■ ■ ■ dxN ipN,t{x, X2,..., XN)i’N,t{y^ X2,..., xn)- 

Notice that we chose the normalization Tr 7 jJ^^ = N. 
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Analogously, for k = 2,3,... , N, we can define the fc-particle reduced density associated with 
'ipN,t by 

The integral kernel of the /c-particle density matrix is given by 

f - (2.5) 

I dxfc-|_i... dx]\f ipNjtixi, ■ ■ ■, Xk, 3Jfc-(-i, • • •, xjsf^Tp /v", 4(2/1 ) ■ ■ ■) 2 /fc) ,..., Xjv) ■ 

The normalization is such that Tr . 

(k) 

Clearly, for k < N, the /c-particle reduced density does not contain the full information 
(k) 

about the system. Still, 7 )^ ^ is enough to compute the expectation of any /c-particle observable: Let 

be an operator on the one-particle space and denote by = 1 (g) • • • (g) (g) • ■ ■ (g) 1 the 

operator on acting like j7) on the i-th particle and trivially on the other {N — 1) particles. 

We write Then 

N 

(V’4v,4,dr(j(b)V;^_,) = = NTrY'^\^N,t){M = Tr 

i=l 

Similarly, if is an operator on the /c-particle space and if we denote by the 

operator acting like on the k particles ii,... ,ik, we have 

{^1 } 

where the sum on the 1 . h. s. runs over all sets of k different indices chosen among 

{!,...,iV}. 

Convergence of reduced densities. It turns out that reduced densities provide the appropriate 
language to describe the convergence of the many-body quantum evolution towards the Hartree 
dynamics. Under appropriate assumptions on the external potential I4xt and, more importantly, on 
the interaction potential V, one can show the convergence of the reduced densities associated with 
the solution of the Schrodinger equation (j2.2p towards orthogonal projections onto products 
of solutions of the Hartree equation (|2.4I) . More precisely, consider a sequence of wave functions 
G Lg(R^^) with reduced density 7 ^^ satisfying 

^ {N^oo) ( 2 . 6 ) 

for a ip £ L^(M^). Let 'ipN,t = be the solution of the Schrodinger equation (|2.2I) with initial 

data ifN- Then we expect, and under appropriate assumptions on Uextj U we can show that 

(N^oo). (2.7) 

Here ipt denotes the solution of the Hartree equation (|2.4p with initial data ipo = ip. The convergence 
in ( 1 ^ can be understood in the trace-class topology. In fact, since the limit is a rank-one projection. 
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weak convergence implies convergence in the trace norrrO- Moreover, notice that convergence of the 
one-particle reduced density towards a rank-one orthogonal projection also implies convergence of 
higher order reduced densities in the limit N ^ oo (the argument is outlined in |48] . after Theorem 
1 ), 

( 2 . 8 ) 

The convergence here (and in (I2.7ld is for fixed t and k. Eq. (12.71) and (|2.8I) explain in which sense 
one should understand the approximate factorization (I2.3p . (For any operator A we use the notation 
A for its /c-fold tensor product.) 

Some results. The first rigorous results of the form (|2.7I) have been obtained for smooth inter¬ 
actions by Hepp [ID] and for singular potentials by Ginibre and Velo [3S]. For bounded potential, 
a proof of the convergence (EZD, based on the analysis of the so-called BBGKY hierarchy, has been 
given by Spohn [62] (see next paragraph). The BBGKY technique has been later extended to poten¬ 
tials with Coulomb singularity in [301 [9] and, for semi-relativistic bosons, in [26]. A new approach 
giving a precise estimate on the rate of the convergence towards the Hartree dynamics has been 
developed in [56] for potentials with Coulomb singularities (and further improved in [21ji: this ap¬ 
proach, based on the ideas of [IHl [35], will be presented in Section [3l Other bounds on the rate of 
convergence towards the Hartree evolution have been established in [H] (based on ideas proposed in 
|54j f and recently in [3]. In [331132] the convergence of the many-body evolution has been interpreted 
as a Egorov-type theorem. In the authors study the propagation of the Wigner measure in 

the bosonic mean-field limit. Next order corrections to the Hartree dynamics have been considered 
in [351 [DD], leading to a better approximation of the many-body evolution. (A related problem is 
the study of the fluctuations around the Hartree evolution, which will be discussed in Section [31) 
Instead of a fixed interaction V, it is also interesting to consider Y-dependent potentials, scaling like 
Vn{x) = Y3“K(Y“ x) (in the three dimensional case) and converging towards a delta-function in the 
limit of large N. In this case (assuming a < 1; for a = 1, one recovers instead the Gross-Pitaevskii 
regime, which will also be discussed in Section [3]), the many-body evolution can be approximated by 
a nonlinear Schrodinger equation with a local cubic nonlinearity. Results in this direction have been 
obtained in HI El ED] in the one-dimensional setting, in [ 32 [ in the two-dimensional case and in m in 
three dimensions. It is also possible to start from Hamiltonians with three-body interactions; in this 
case, the evolution can be approximated by a quintic nonlinear Schrodinger equation; see [331 [T5] . 


The BBGKY approach. The main idea of the BBGKY approach, which was first applied to 
many-body quantum systems in the mean-field regime in [32], is to study directly the evolution of 
the reduced densities defined in (|2.5p . To explain this idea it is convenient to normalize the reduced 
densities associated with the solution of the Schrodinger equation, defining, for /c = 1,..., Y, 



1 (fc) 
'N\ '^N,f 
k) 


The new density matrices are normalized so that Tryj^') = 1 for all Y € N and all A: = 1,..., Y. From 

the Schrodinger equation for V’iv.t it is easy to derive evolution equations for the family 

It turns out that the evolution of the Y reduced densities is governed by a hierarchj0 of Y coupled 


®First, by testing the difference N ^7^) — \ipt){‘Pt\ against \ipt){'Pt\ it implies Hilbert-Schmidt convergence. Then, 
since \ipt){<pt\ is a rank-one projection, the operator \ has exactly one negative eigenvalue. (If there 

were two linearly independent eigenvectors ^1,^2 with negative eigenvalue, one could find a linear combination ^ such 
that < 0 .) Since Tr7)^j — \ipt){ipt \ = 0 , its absolute value is equal to the sum of all positive eigenvalues, and 

therefore the trace norm is twice the operator norm. 

^i. e. the equation for 7^) depends on 
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equations, known as the BBGKY hierarchy (BBGKY stands for Bogolinbov-Born-Green-Kirkwood- 
Yvon): 


j=l i<j 


+ ^ ^Trfc+1 V{xj - 


N 


(fc+i)' 

t 


J = 1 


(2.9) 


where we nse the convention that = 0 and where Tr^+i denotes the partial trace over the 

degrees of freedom of the {k + l)-st particle. Here we also introduced the commutator, defined by 
[A, B] = AB — BA on a domain suited to the operators A and B. Notice that the second term on 
the r. h. s. of (12.9p describes the interaction among the first k particles in the system. The last term 
on the r. h. s. of (12.9p . on the other hand, corresponds to the interaction of these k particles with 
the other (N — k) particles. At least formally, the BBGKY hierarchy (|2.9I) converges, in the limit of 
large N, towards the infinite hierarchy 




K 

Y[ 

i=i 


-A, 


~{k) 

i ’ 7oo,t 


+ ^Trfc+i V(xj - Xk+i),7, 


i=i 


(k+i) 

oo,t 


( 2 . 10 ) 


It is simple to check that this infinite hierarchy has a factorized solutions given by 

products of the solution of the Hartree equation (|2.4p . This observation suggests a general strategy 
to show the convergence (12.8p of the reduced densities towards projections onto products of solutions 
of the Hartree equation. The strategy consists of three steps: 

• Compactness: first, one needs to prove the compactness of the sequence (in N) of families 

Tw,* = with respect to an appropriate weak topology. Gompactness implies in 

particular the existence of at least one limit point Too,* = 

• Convergence. Secondly, one needs to characterize limit points of the sequence Tas solutions 
of the infinite hierarchy (I2.10p . In other words, one has to show that any limit point Too^t of 
the sequence T^v,* satisfies ()2.10p . 

• Uniqueness. Finally, one has to prove the uniqueness of the solution of the infinite hierarchy. 
This implies immediately that the sequence T 7 v,t converges, since every compact sequence with 
at most one limit point converges. Moreover, since we know that the factorized densities 

= \ipt){(pt\®^ are a solution of (|2.1Up . uniqueness also implies that 7N,t for 

all A: G N (the argument proves convergence with respect to the weak topology with respect 
to which one showed compactness in the first step; however, since the limit is a rank-one 
projection, weak convergence immediately implies convergence in the trace norm). 

The most difficult of the three steps is the proof of the uniqueness of the solution of the infinite 
hierarchy. Let us illustrate how to prove uniqueness in the case of a bounded interaction potential 
V G L°°(M^). To this end, we rewrite the infinite hierarchy (|2.I0p in integral form as 

7S=«*'‘’(‘)77o + - (2.11) 

^ JO 
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Here we defined the action of the free evolution (t) on a A:-particle density 7 ^*^^ by 

and the action of the collision operator on a {k + l)-particle density 7 !^+^) by 


Bik)^{k+i) = ^Trfc+i -xfc+i),7("+') 

i=i 


( 2 . 12 ) 


Notice that the collision operator maps {k + l)-particle density matrices into /c-particle density 
matrices. It is important to observe how the free evolution and the collision operator affect the trace 
norm. On the one hand, we clearly have 


On the other hand, for bounded interaction potentials, we find 

k 

^ [tv \v{xj - Xk+ih^’^+^^\+TT\j^'^+^'^V{xj - Xk+i) 
i=i 

<2fc||H||ooTr|7('=+i)| 

<2fc||H||oo||7^"+'^l|tr. 


(2.13) 


(2.14) 


Here we used that Tr|Hil| < ||H||TV|H| for any bounded operator A and any trace-class operator B. 
(In the same way Tr|HH| < ||i?||Tr|H| if B is bounded and A trace class.) 

Iterating (12.111) . we obtain the n-th order Dyson series 


'OO,0 

-\ rt r 

+ E^/*—/ 

* Jo Jo 


^m — 1 


dSr 


m=l 


X {t - Si)H W • • • U^k+m-l) 


+ 


- f 

Jo 


^^n—1 


dsi ... dSn 

Jo 


Suppose now that the families ri^oo,t = { 7 ! t}k>o and r 2 ,oo,t = {72 cL t}k>o are two solutions of the 
infinite hierarchy, written in the integral form (12.IIP , with the same initial data. Expanding both 
solutions as in (I2.15p . and taking the difference (so that all fully expanded terms cancel), we find 


Tr 


~(k) ~{k) 

7l,oo,t 72,oo,t 


/*t pSn — i 

Jo Jo 

pt pSji—\ 

+ *„ ||M<‘>(t - si)Bl‘)... B<‘+"-i)7<‘+2 

Jo Jo 

Applying iteratively the bounds ()2.13p and (I2.14p . we conclude that 
Tr 


|tr- 


(fc) (fc) 

Tl,oo,t 72,00,4 


< 2^k{k + l)...{k + n- l)(2||H||oo)" < 2"(4||H||oo|t|)", 
n! 
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where we used the normalization Tr 


(k+n) 


rr« (fc+n) 

72,00,/ 


1 and the bound 


fn + k -l\ 

V k-i J-^ 


For |t| < (8||F||oo) we obtain that for all n G N 


Tr 


~(k) 

'yi,oo,t 


~(k) 

■ 7w 


< 2 


k—n 


^(k) ~(’fc) 

Since the 1. h.s. is independent of n, it must vanish. This proves that 7i = 72 <L t rl — 

(8||F||oo)~^- Since this argument only uses the normalization Tryj^^”^ = Try^^^”^ = 1, which hold 
for all t G M, it can be iterated to prove uniqueness of the solution of the infinite hierarchy for all 
t G M. 

For a Coulomb potential V{x) = ±l/|a:|, the proof we outlined above can be modified by intro¬ 
ducing a different norm for density matrices (this approach was first used in [30]). For a /c-particle 
density 7 ^^/ acting on we define the Sobolev-type norm 




= Tr 


5i... ...Si 


(2.16) 


where S'j = (1 — Since the Coulomb potential is bounded with respect to the kinetic energy 

in the sense that as operators 

±rT < Cil - A) 

|x| 

one can show that the collision operator defined as in (I2.12p . but now with V{x) = ±l/|x|, 

C Q "i"! ch OC 

This bound replaces (12.1411 . Using this new bound one can prove, similarly as explained above for 

the case of bounded potentials, that any two solutions Fi^oo^t = t\k>i and r 2 ,cxD,t = {n^lo t\k>i 

of the infinite hierarchy with the same initial data satisfy 


~(fc) ~(fc) 

7i,oo,t - 72,oo,t 


H- 


(k) - 


< c 


k—n 


^(fc+ri.)|| ||~(fc+ri.)|| 

yi.^oo,t 11 ^ 2 , 00 ,^ 


To conclude uniqueness, here we need to show a-priori estimates of the form 


~{k) 

7 oo[t 


H 


(k) - 


< C'' 


(2.17) 


for all t G M and /c G N, valid for any limit point Foo,* = {7 (^\}a:>i of the sequence of families 

F N,t = { 7 ^!}^! of densities associated with the solution of the Schrodinger equation. (This step 
was not needed in the case of bounded potentials because the trace norm trivially remains uniformly 
bounded.) To obtain a-priori bounds of the form ()2.17p . one can use energy conservation; for details, 
see [30] . 
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3 Coherent States Approach 


In this section, we illustrate the method developed in [SB], based on the original ideas of [IBIES], to 
prove the convergence of the many-body evolution towards the Hartree dynamics. This method is 
based on a representation of the many-boson system on the Fock space and on the study of the time- 
evolution of a special class of initial data, known as coherent states. With respect to the BBGKY 
approach that was presented in the last section, the analysis of the evolution of coherent states allows 
us to obtain precise bounds on the rate of convergence towards the Hartree dynamics. This is an 
important point, since real systems have a large but, of course, finite number of particles. Bounds 
on the rate of convergence are therefore crucial to establish whether the Hartree equation is a good 
approximation to the dynamics of a given boson gas. 

Fock space. The bosonic Fock space over is defined as the direct sum 

.F=0L2(M3-), 

n>0 

where denotes the subspace of consisting of all permutation symmetric functions. 

On T", we can dehne the inner product 

n>0 


We denote by II.Hj- the corresponding norm. We always use vectors T G T" normalized such that 

n>0 


The Fock space allows us to describe states of the system where the number of particles is not 
fixed. The vector T = £ J- describes a state with probability ^ particles. In 

particular, states with exactly N particles are embedded in the Fock space; they are described by 
vectors {0,0,..., ipNjO, ■ ■ ■} £ having only one non-zero component. These vectors are eigenvectors 
of the number of particles operator M defined by 


for any T = such that 

< oo. 

n>0 

The vacuum vector n = {I,0,...} gJ-” plays a special role; it is an eigenvector of Af with eigenvalue 
zero and describes a state with no particles at all. 

Creation and annihilation operators. For any one-particle wave-function / G we define 

the creation operator «*(/) and the annihilation operator a{f) by setting 


1 

(a*(/)T)(’^)(xi, . . . , Xn) = ^ V (xi, . . . , ...,Xn) 

,=1 

(o(/)^')(”)(xi, . . . ,Xn) = Vn + 1 J dxf{x)'tlj^'^~^^^{x,xi,...,Xn). 

The interpretation is straightforward: a*{f) creates a new particle with wave function /, while a{f) 
annihilates such a particle. Creation and annihilation operators are closed densily defined operators 
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on T\ moreover, a*(/) is the adjoint of a(/) (as the notation suggests). Notice also that a*(/) is linear 
in /, while a(/) is antilinear. Creation and annihilation operators satisfy canonical commutation 
relations 

[a{f),a*{9)] = {f,9), [«(/), 0 ( 5 )] = [a*{f),a*{9)] = 0. (3.1) 

Despite the bosonic creation and annihilation operators being unbounded, usually domain questions 
are unproblematic since one can take the domain of a sufficiently large power of the number operator 
to easily make sense of most expressions and in particular of their commutators — see (|3.7I) . (One has 
to be more careful in this respect when working with the operator-valued distributions below.) It is 
also useful to introduce operator-valued distributions a* and which formally create and annihilate 
a particle at position x G M^. It terms of these distributions 


a*{f) = J f{x)aldx, a{f) = J f{x) a^dx 


and 


— d{x y), [ua;,Q-j/] — [a 2 ,,aj^] — 0. 

Notice that, if 'k = 

{ax^)^"‘'>{xi, ...,Xn) = VnTl'ip^'"'^^\x,xi ,... ,x„). 

Hence 

(^', AA$) = ^ n f dxi... dxncp^"''^ (a^i, • • •, (xi,..., x^) 

n>0 

= ^ / dxdxi . . . dXn-l{ax^f^ ^\xi, . . . ,Xn-l){ax^)^'"~^\xi, . . . ,Xn-l) 

n>l 

= / dx{ax^,ax^) 


(3.2) 

(3.3) 


and in this sense we can write 

Af = J dxa*ax- (3.4) 

This expression for Af suggests that, although creation and annihilation operators are unbounded 
operators, they can be bounded with respect to the square root of the number of particles operator. 
We have 

±(T,(a(/) + a*(/))T) <2KT,a(/)T)| 

< 2 I dxdxi... dxn'ip^'^'^ {xi ,..., Xn) f{x) (x, Xi, . . . , Xn) 

n>0 

<2J;^^^TT||/||||V^(")||||V^("+')|| 

n>0 

(3.5) 

and therefore in the sense of forms 

±(a(/) + a*(/))<4||/||Ari/2. ( 3 . 6 ) 
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In norm, we have the bounds 


ll»(/)'I'lb < ll/lll|V‘''^4'|b and 

for any / G L^(R^). To prove the first bound in (|3.7I1 . we observe that by the Cauchy-Schwarz 
inequality 

^< 11 / 11(1 cixiia^n^)' m 

The second estimate in (13.7p follows from 

||a*(/)T||3. = (T,a(/)a*(/)T) = (T, a*(/)a(/)T) + \\ffm& = ll«(/)'I^f + ll/f ll^f 

and (|3.8I) . 

Similarly to (13.4p . we can express the second quantization of any one-particle operator in terms 
of the operator-valued distributions a* and a^: Let be an operator on the one-particle space 
The second quantization of is the operator on T" defined by the requirement 

that 

n 

E 

i=l 

where denotes the operator acting on as J(i) on the Tth particle and as the identity 

on the other (n — 1) particles. If the one-particle operator has the integral kernel y), we 

can write 


|a(/)^||j- < j dx\f{a 




n>l j=l 

= ^n/ dxdydx2 ■ ■ ■dxnip^"'Hx,X2, ■ ■ ■ ,Xn)j'^^Hx;y)'tp^‘^'>{y,X2,... ,Xn) 

n>l d 

= ^ / dxdydx2 . ..dxn y)(aa; 4 >)("--L(x 2 , ...,Xn) (oyT)(”'“^)(x2,... ,x„) 

n>l d 


J dxdy J^^\x;y) {ax^,ay^). 


Thus 


dT{J^^^) = J dxdy J^^\x\y)L 


(3.9) 


Since the number of particles operator is the second quantization of the identity, N = dr(l), the 
last expression is consistent with (13.4p . Notice that if is bounded, then its second quantization 
is (although generally unbounded) bounded with respect to the number of particles operator, 

i.e. 


|(V^,dr(j(^))V^)| < ||j(^)||(V’,AAV’) for all V’ € .F. (3.10) 


and also 

||dr(j(i))T|| < ||jW||||A^T|| . 
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We can extend (|3.9I) to operators involving more particles. Let be an operator on the 
A;-particle Hilbert space We define the second quantization ) of by 

(dr(jW)^)W= (3-11) 

where the sum runs over all sets {zi,..., of k different indices in { 1 ,..., n}, and where 
denotes the operator on acting as on the particles ii,... ,ik and as the identity on the 

other (n — k) particles. If has the integral kernel ... ,Xk;yi, ■ ■ ■, yk), we find, similarly 

to (|3.9p . 

dr{J^^'^) = J dxi...dxkdyi...dykJ^'"\xi,...,Xk-,yi,...,yk)al^ .. .al^Uy^.. .ay^. (3.12) 

Similarly to (I3.10p . we find 

|(V’,dr(j(^))V^)| < ||j(^)||Ar(Ar-l)...(AA-A; + l). 

Reduced densities. We define the one-particle reduced density associated with a normalized Fock 
space vector 'I' G as the operator 7 ^^^: L^(]R^) —>• L^(]R^) such that 

Tr j(i)^(i) = 

for all one-particle observables From (13.9h . we find that the integral kernel of the one-particle 
reduced density 7 ^^^ is given by 

7 (^)(x; y) = (3.13) 

We define the fe-particle reduced density associated with the Fock space vector 'k G as the operator 
7 ^^) on L^(M^^) such that 

Tr = {^,dr{J^’^^)^) 

for every /c-particle observable . Using (|3.12l) we hnd that the integral kernel of the /c-particle 
reduced density 7 ^^^ is 

j’'’^\xi,...,xk;yi,...,yk) = 

Notice that 7 ^*^^ is normalized such that Tr 7 *^^) = ('I',A^(A^— 1)... {J\f — k + I)'!’). 

Representation of the Hamiltonian in second quantization. We define the Hamilton operator Hm 
on the Fock space J-” by with 

n 1 

= ^-^xi + j^^V{xi-Xj). (3.14) 

j=l i<j 

Since it does not play an important role in our analysis, we neglect here the external potential (but 
we could also add it, requiring only straightforward changes). The operator "Hat leaves the number 
of particles invariant, i. e. it commutes with the number of particles operator M. When restricted to 
the Wparticle subspace (consisting of the vectors of the form { 0 , 0 ,..., 0 , 'tpNjO, ■■■})) coincides 
precisely with the Hamilton operator (|2.ip defined in the previous section (of course up to the external 
potential that we neglect here). Notice that the parameter N appearing in the Hamiltonian (I3.14p 
is not yet related with the number of particles in the system which in Fock space is arbitrary. Of 
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course, later on we will restrict our attention to initial data in T that has approximately N particles; 
this is important to make sure that we are in the same mean-field regime discussed in Section [21 
We can write the Hamiltonian in terms of the operator-valued distributions and a*. From 

(|3.3p we find that 

T-Ln = J + ^ j dxdyV {x - y)ala*yayax ■ (3.15) 

The r. h. s. should be understood in the sense of forms, similarly to ()3.4p . The fact that %m commutes 
with the number of particles operator N is evident from (I3.15p . because, in every summand, the 
number of creation operators matches exactly the number of annihilation operators. In particular, 
it follows that 

IN) 

..., 0, V'iv, 0,... } = {0,..., 0, Q,... }. 

On the r. h. s. we recover exactly the mean-field evolution discussed in the previous section. So what 
have we gained by switching to the Fock space representation of the bosonic system? The answer 
is that now we have more freedom in the choice of the initial data; in particular, we can consider 
initial data having a number of particles which is not fixed, i. e. is a superposition of different particle 
numbers. We will make use of this freedom by considering a special class of initial data known as 
coherent states. 

Coherent states. For / G we introduce the Weyl operator 

W{f) = 

The coherent state with wave function / G is defined a^ 

W{m = ^ g-||/||V2 |i^ ^ ^ I . (3.16) 

Since W{f) is unitary, 

w{fr = w-\f) = w{-f), 

coherent states are always normalized, ||kF(/)n|| = 1 for all / G L^(M^). Weyl operators generate 
shifts of creation and annihilation operators, i. e. 


Wifra{g)W{f) = a{g) + {gJ), 
W{fra*{g)W{f) = a*{g) + {f,g). 


This implies that coherent states are eigenvectors of all annihilation operators, i.e. 


a{g)W{f)n = W{f)W*{f)aig)W{m = W{f) {a{g) + {gj))n = {f,g)W{f)n. 

From (|3.16D . we see that coherent states do not have a fixed number of particles. Instead they are a 
linear combination of states with all possible numbers of particles. The expectation of the number 
of particles operator in the state W{f)Q is 

= J dx{wim,aia,w{m) 

= J dx{W{m,{al+7ix)){a, + fix))n) = \\ff. 

®To evaluate W (/)n we made use of the Baker-Campbell-Hausdorff formula for operators A, B with the property 
that [\A,B],A] = [[A,B],B] = 0; 

^A+B 

Since a{f)Q = 0 one then simply has to expand the exponential of a*{f). 
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More precisely the number of particles in the coherent state W{f)Q is a Poisson random variable 
with mean and variance ||/|p. 

Dynamics of coherent states. We are interested in the evolution of coherent initial states. For 
(f E with lly^ll = 1, we consider the coherent state W{s/N ip)D. The expectation of the number 

of particles operator is N] for this reason we expect to be in the mean-field regime, in which many- 
body interactions can be effectively approximated by an average potential. This is made rigorous by 
the following theorem. 

Theorem 3.1. Let V be a measurable function satisfying the operator inequality 

p2(a;) < C(1 - A) (3.18) 

for some C > 0, and let ip E Let {.^ivjArgN be a sequence of vectors in F with ||.^iv|| = 1 and 

such that for some C > 0 

< c VA E N. (3.19) 

Let be the one-particle reduced density associated with Then there 

exist constants D,K > 0 such that 


Tr 






(3.20) 


for a//1 E M and all N gN. Here pt denotes the solution of the nonlinear Hartree equation 

idtpt = -A<y9t + (P * \ptf)pt (3.21) 


with initial data po = p. 

Remarks: 

• The trace-norm bound ()3.20p is of order one; this is to be compared to TryjJ^^ = TrN\pt){pt\ = 
N. In this sense, (|3.2UI) establishes a relative rate of convergence of order 1/N. 

• The assumption V‘^{x) < C{1 — A) holds for potentials with a Coulomb singularity, since 
Hardy’s inequality implies that 

In fact, the result (j3.20l) can be easily extended to potentials V E L^(]R^) -|- L°°(]R^), although 
in general the time dependence on the r. h. s. of (13.2011 may be worse. 

• The condition (I3.19P guarantees that the initial deviations from the coherent state W{y/Np)Ll, 
which are described by the vector are small compared to the number of particles in the 
initial data W{^/Np)^N, which has expectation N. 

Sketch of the proof of Theorem \3.1\ . According to (I3.13P , the one-particle reduced density associ¬ 
ated with the Fock space vector Thas the integral kernel 

y) = a*yax'^N,t)- (3.22) 
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We define the fluctuation vector at time t by setting 

^N,t = e-^^^^W{^/Nip)CN = W{VNipt)CN,t (3.23) 

where ipt is the solution of the Hartree equation (|3.21l) . Equivalently, we can write 

iN,t = UN{t] 

where we introduced the fluctuation dynamics 

UN{t-, s) = W{y/Nipt)*e-^'^^^W{y/N^s)- (3.24) 

Plugging the r. h. s. of (I3.23P in (|3.22l) . we obtain 

lN]ti^-,v) = {W{VNn)^N,t,ala^W{VN^t)^N,t) 

= {a*y + '/Nlpt{y)){a^ + '/N(pt{x))^N,t) 

= N'^t{y)n{x) + '/N^t{y){CN,t,ax^N,t) + VN(ptix){CN,t,alCN,t) + (Cat,*, OyOx^Af,*). 

Integrating against a one-particle observable with kernel J^^\x]y), we find 

= VN(^^N,t, [a(j(^Vt) + a*(^^'Vt)] ^N,t) + (3.25) 

According to (13.10^ and (13.611 . the expectations of the operators a(J^^Vt) +o*(J^^Vt) and dr(j(^)) 
can both be bounded by the expectation of the number of particles operator, in the state = 
U]y(t‘,0)^N- The next proposition is taken from [T71 Prop. 3.1] but similar estimates have already 
been proven in [561121] . The proposition extends also to higher moments of the number of particles 
operator. 

Proposition 3.2. Let V be a measurable function satisfying the operator inequality \3.18\) and let 
ip E i/^(M^). Then there exist constants D,K > 0 such that 


{lLNit-,0)'il>,MUNit-,0)'il>) < {f), 

for any fj £ T. Here Ujy is defined as in ^3.24\)- 


N+-N^ 


f}) 


While Proposition 13.21 immediately implies the desired bound for the second term on the r. h. s. 
of (13.251) . some more work is required to get rid of the additional y/N factor in the first term on the 
r. h.s. of (|3.25p . We will omit the details which can be found in |56[ I21j . and just write the final 

1 m cs Eo 

Since the space of trace class operators /1^(-L^(M^)), equipped with the trace norm, is the dual of the 
space of compact operators, equipped with the operator norm, the last bound implies (|3.2ni) . 

Sketch of the proof of Provosition \3.2X Let us briefly discuss the main ideas needed to bound the 
growth of the number of fluctuations. We notice, first of all, that the fluctuation dynamics Z^ 7 v(f; s) 
satisfies the Schrodinger type equation 


idtUN{t] s) = CN{t)UN{L s) 
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with the time-dependent generator 


CN{t) = \idtW*{y/Nipt) \ W{y/Nift) + W*{y/Nipt)'HNW{y/Nipt)- 
The time derivative i^ 

idtW*{^/N^t)\ W{^/N^t) = Ci^N{t) -^/N{a*{idm) + a{idm)) 

and, using (I3.17jl . that 

W*{^fN^t)nNW{^^t) = C2,N{t) + y/N {a*{-Aipt + * \^t\^)ipt) + a{-A^t + * Wt?)^t)) 

+ j dxVxa*Vxax + j dx{V *\ipt\'^){x)alax 

J dxdyV(x - y)Lpt{x)Tp^{y)alay 

J dxdyV(x - y) {(ptix)ipt{y)alal + lpi-{x)i^t{y)a^ay) 

j dxdyV(x - y)a* {iptivWy + ^tivhy) 

^ J dxdyV{x-y)c 


+ 

+ 

+ 

+ 




for appropriate constant^ C'i^Ar(t) and C 2 ,N{t)- Since satisfies the Hartree equation, the contri¬ 
butions of order y/N (which are linear in a* and Ox) cancel exactly, and therefore (up to a constant 
contribution, which can be absorbed in the phase) we obtain 


^Nit) = j dxVxalVxax + j dx{V * \ipt\‘^){x)alax + j dxdyV(x - y)(ptix)(pt{y)a*ay 


+ J dxdyV{x-y) {(pt{x)^pt{y)ala*y + ^p^{x)^p^{y)axay) 
1 


+ 


y/N 


J dxdyV{x-y)al ((^t(y)a* + (^t(y)ay) 


(3.26) 


1 


dxdyV{x - y)ala*yayax. 


We notice that C]yf{t) contains terms (the second and third line) which do not commute with the 
number of particles operator. This means that, as expected, the fluctuation dynamics 7/Ar(t), in 

®A systematic way to calculate to calculate time derivatives of the form (with A(t) a sufficiently 

regular family of operators) is as follows. Start by writing 

= lim i = lim i [ dX-^ (^-A(t+h)\^A(t)\\ 

h^o h \ / h^o h Jq dX \ / 

= -lim r dX ~ dA 

^ Jo 

In the integral then use the Baker-Campbell-Hausdorff formula e~^Be^ ~ ^ ~ Jo dpe~^'^[A, B]e^^ for operators A, B. 
The commutator [A, B] in our case is a complex number, so the expression is easy to evaluate. (Otherwise it could be 
helpful to iterate the last integral formula to obtain a series.) 

^We speak of constants since these are not operators, but rather just complex numbers. Of course they depend on 
time and on N. 


19 









contrast with the original dynamics e does not preserve the number of particles. Nevertheless 

we can control the growth of the expectation of Af. We compute 

idtiUisfit; 0)^N0 )Cn) 

= [UNit] 0)^Ar, [M, CN{t)\UN{t-, 0),^Ar) 

= 2ilm j dxdyV{x - y)(ptix)(pt{y){UN{t-,0)^N,[Af,a*a*y]UNit;0)^N) 

J dxdyV{x - y)ipt{x){UN{t;0)^N, [Af, a*a*yax]UNit]0)^N)- 

Using the canonical commutation relations (13.21) we find 


2i ^ 

H-;=lm 

Vn 


dt{ldN{t;0)^N,AfUN{t-,0)^N) = 4Im 


dxdyV{x - y)ipt{x)(pt{y){UN{t;0)CN,al.a*yUN{t;0)CN) 


H—;=Im 

Vn 


J dxdyV{x - y)ift{y){ldN{t]0)^N,ala*yaxUN{t;0)^N) 


=: I + II. 


(3.27) 


The first term can be bounded using (|3.7I) : 


HI = 4 


< 


< 


J dxipt{x){axUN{t]0)^N,a*(y{x - ■)(pt)UN{t]0)^N) 

4 J dx\ipt{x)\\\ao,UN{t;0)^N\\\\a*{V{x - ■)ipt)UN{t-,0)^N\\ 

4 J dx\ipt{x)\\\axUN{t;0)^N\\\\V{x - ■)ipt\\2\\{Af + iy^‘^UN{t;0)^N\\ 


< 4sup ||U(x - ■)tpt\\ 2 \\{Af + l)^/^ZYAr(t;0)^Ar|p. 

X 

From the assumption (14.41) . we obtain 


\\V{x - ■)iptf = / dyV'^ix-y)\ipt{y)\‘^ <C\\ipt\\Hi. 


Furthermore, the Hartree energy 

£H{y^) = j \V<p{x)\‘^dx + ^ j dxdyV{x - y)\(f{x)\‘^\(f{y)\‘^ (3.28) 

is conserved along the Hartree evolution, i. e. for a solution of the Hartree equation (pt, is 

independent of t. Together with the assumption (14.41) . this implies that there exists a universal 
constant C > 0 such that 

llT’dl/fi < C'IIt’oIIj/i- 

Hence 

sup ||U(x - ■)pt\\ < K 

X 

for a constant K > 0 depending only on the constant CO in ()4.4p and on the H^-norm of the initial 
wave function pQ. We conclude that 

|I| < 0)^N, {Af + l)Z^Ar(f; 0)(n)- (3.29) 
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To control the second term on the r. h. s. of (|3.27l) . we proceed as follows: 
2 


|n| < 


< 


Vn 

2 

y/N 

K 


J dx\\a{V{x - ■)(ft)axl{Nit]0)(N\\\\axldNit]0)(N\\ 

sup ||y(x - •)y?i ||2 / dx\\axAf^^^HNit]0)^N\\\\axUNit;0)^N\\ 

X J 




In the first term on the r.h. s. of the last equation, we are going to use the 1/N factor to reduce 
the exponent of the number of particles operator. Since U^{t; 0)MUNit] 0) measures the number of 
fluctuations, it is heuristically clear that it can be bounded by the total number of particles N] more 
precisely, one can prove the operator inequality (see m, proof of Prop. 3.1) 

1 Af‘^ 

-uUt;OW^UNit-,o)<UMt;0)J\fUNit-,o) + 

Applying this bound for we conclude that 

— (Z^Ar(t;0)^Ar,AA^Z^Ar(t;0)^Ar) < (Z^Ar(t; 0)^Ar, 0)^Ar) + — (^at, AA^^at) 

and therefore that 

|II| < K{UN{t-, 0)Civ, (AA + l)UN{t-, 0 )^n) ■ 

Inserting the last bound and (|3.29D into (|3.27p . we find 


— 0)^n) 


< K{UN{t; 0)^N, (AA + l)UN{t; 0 )^n)- 


Gronwall’s Lemma implies therefore 


This concludes the sketch of the proof of Proposition 13.21 more details can be found in [MIEI]. 


4 Fluctuations around Hartree dynamics 


The coherent state approach presented in the last section also allows us to describe the fluctuations 
around the Hartree dynamics. 

Fluctuations for coherent initial states. For simplicity, let us consider the evolution of a coherent 
state initial data Tat = W{y/Nip)fl E T. (As in the last section, we could also consider approximate 
coherent states Tat = with ^Af satisfying (^at, (AA^/^ + A^“i/2y7^/2 _|_ j 

According to the definition of the fluctuation dynamics (j3.24D . we can write 

e-^'^^*W(VNip)Q = W(VN(pt)ldN(t;0)n. 


Recall the generator CN{t) of the fluctuation dynamics UjyfiFs) as given by (13.261) . The last two 
terms on the r. h. s. of (j3.26p . the cubic and the quartic contribution to C]sf{t), seem to vanish, in the 
limit of large N. Therefore we define a new time-dependent generator 


Coo{t) = j dxVxalVxttxF J dx{V*\(pt\‘^){x)a*ax + J dxdyV{x - y)<ptix)(Pt{y)c 
+ / dxdyV{x-y) {^pt{x)ipt{y)ala*y+Tpt{x)Tp^{y)axay) 


(4.1) 
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keeping only the quadratic part of C]\f{t). We denote by Uooit', s) the evolution generated by Coo{t), 

i.e. 

idtUoo{t\ s) = Coo{t)Uoo{t] s), Uoo{s\ s) = 1. (4.2) 

For interaction potentials satisfying the condition (I3.18|) one can prove that 


\\{UN{t;0) -U^{t;0))iP\\ < 


c\A 

Vn 


||AA3/2V>|| + 


y/N 


IIA^VII 


(4.3) 


for some constant C > 0. The proof of (14.3p is based on the identity 

UN{t;0) -Uoo{t;0) = f UN{t] s) {Cn{s) - Coo{s))Uoo{s;0) ds 
Jo 

and on the control of the growth of moments of the number of particles operator with respect to the 
limiting fluctuation dynamics Woo(s;0). The bound (14.311 leads to the following theorem. 


Theorem 4.1. Let V be a measurable function satisfying the operator inequality 


V^ix) < C(l-A) 


(4.4) 


for some C > 0, and let ip E 


Then there exists a constant D > 0 such that 


e-^'^^^W{y/Np)VL - W{y/N pt)Uoo{t-,t))VL 


< 


m 

y/N' 


In other words, Theorem |TT] states that, if we approximate the many-body evolution of the initial 
coherent state taking into account also the limiting fluctuation dynamics Uoo{t]0) rather than only 
the Hartree dynamics of the coherent state, we get convergence in norm as N ^ oo (with a rate 
proportional to while approximation of the many-body evolution just by evolved coherent 

states only gives convergence in the weaker sense of reduced densities. 

Bogoliubov transformations. In contrast with ILjyf^t] s) the limit fluctuation dynamics Uoo{t', s) is a 
simple evolution, because its generator is quadratic. In fact, it acts as a Bogoliubov transformation, 
a general concept that we introduce now. For f,gG L^(M^), we define 


= a{f) + a*{g). 


Then we have 

A*{f,g) = a* if) + a{g) = A{g,J) = A{J{f,g)) 
where J : L^(M^) 0 L^(R^) —)• L^(M^) 0 L^(M^) is the antilinear map defined by 


J{f,9) = (gJ)- 


The canonical commutation relations take the form 

[A{fi,gi),A*{f2,g2)] = {{fi, gi), S{f2, g2)) 

with the linear operator S : L^(M^) 0 L^(M^) —)• L^(M^) 0 L^(M^) defined by 


S{f,g) = {f,-g), 


(4.5) 


(4.6) 
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and where (•, •) denotes the inner product in A Bogoliubov transformation is a 

linear map 0 ; © L^(M^) —)• © L^(M^) that preserves (|4.5p and (|4.6I) . In other words, 

it is a linear map satisfying the two conditions 

QJ = JQ and 0*S0 = S. 

It is easy to check that 0 is a Bogoliubov transformation if and only if it can be written in the block 
form 

--['il) (-) 

with two linear operators U, V : L^(R^) L^(M^) satisfying U*U — V*V = 1 and U*V + V*U = 0. 

A Bogoliubov transformation 0 : © L^(M^) —>• L^(R^) © L^(R^) is called implementable if 

there exists a unitary map ^ : J- ^ J- defined on the bosonic Fock space J- over L^(R^) such that 

E*A{f,g)E = A{e{f,g)) 

for every f,g^ L^(R^). The Shale-Stinespring condition states that a Bogoliubov transformation 0 
is implementable if and only if the off-diagonal operator V : L^(R^) —)• L^(R^) is Hilbert-Schmidt 
(i.e. if TtV*V < oo). 

The limit fluctuation dynamics hloo{t;s), defined in (I4.2p . acts for any t,s G R as a Bogoli- 
ubov transformation: It is simple to check that there exists a two-parameter group of Bogoliubov 
transformations 0(t; s) : L^(R^) © L^(R^) —>• L^(R^) © L^(R^) such that 

s)A{f, g)U^{t-, s) = A(0(f; s){f, g)). (4.8) 

The time-dependent Bogoliubov transformations Q{t;s) are easily seen to satisfy the equation 

idtQ{t\ s) = D{t)Q{t\ s) 


with the generator 

—A + (I^ * + Ai A 2 

A 2 —A + (I^ * + Ai 

Here we denoted by Ai and A 2 the operators with integral kernels Ai(x;y) = V{x — y)(pt{x)ipt{y) and 
A 2 {x; y) = V{x — y)ipt{x)^pt{y) (these are of the form of the exchange operator, but their physical role 
is different). Identifying Uoa{t; s) with a Bogoliubov transformation means that the limit fluctuation 
dynamics can be determined by solving a partial differential equation on L^(R^) © L^(R^); in this 
sense we can think of Woo(t; s) as a simple effective dynamics (in contrast with s) which a true 
many-body evolution). 

Probabilistic interpretation. The convergence of the reduced densities can be interpreted, in the 
language of probability theory, as a law of large numbers. Let be a one-particle observable, i. e. 
a self-adjoint operator over L^(R^). We denote by the operator on L^(R^'^) acting as on the 
f-th particle and as the identity on the other (N — 1) particles. W. r. t. a factorized Al-particle wave 
function , the observables define independent and identically distributed random variables, 
and therefore, for every d > 0, 









(99, 



^ 0 


{N —>■ 00 ). 
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W. r.t. the evolved wave function solution of the mean-held Schrodinger equation 


N 


N 


j=l i<j 


with initial data V'iv.o = the N observables • • •, are no longer independent. Neverthe¬ 
less, it is easy to check that the convergence of the reduced densities implies that 








^ 0 


{N oo). 


At time t = 0 we also have a central limit theorem stating that the huctuations of , 

appropriately normalized, are Gaussian in the limit. Does the same hold true for t ^ 0? The answer, 
obtained in mm, is positive; with respect to the measure induced by V'A'",*) we have 


- ((^t, ^ Gauss(0, af) 


(4.9) 


as N ^ oo, in distribution. The variance of the limiting Gaussian variable is given by 


= 


(0(t;O) {Oipt,JOipt),Q{t]0) {Ocpt,JOcpt)) - 


0 (t; 0 ) {Ocpt, JOift ), (<y?, ¥?) 


where 0(t; s) denotes the Bogoliubov transformation dehned in (14.811 to describe the limit fluctuation 
dynamics. Eq. (gJl) shows that in the mean-field limit the correlations among the particles are weak 
enough for the central limit theorem to hold true, but they are strong enough to change the variance 
of the limiting Gaussian variable. (For a completely factorized wave function the variance 

would clearly be ((^i, - ((^i, 

Fluctuations for data with fixed number of particles. For A’-particle initial data, a better ap¬ 
proach to study fluctuations around the Hartree dynamics has been proposed in [46] , based on ideas 
developed in m to analyze the excitation spectrum of mean-field Hamiltonians. Fix ip E L^(M^). 
Then we can write any A^-particle wave function ipN ^ a^ 

= '00) -H '0(1) (g)^ -^ 0(^“l) (g)^ 99 -h 0(^) 

where the n-particle bosonic wave function 0 (’^) is assumed to be orthogonal to 99 in all its n entries. 
This defines a map 

: L2(]R3iV) ^ 

ipN ^ , • • •, , 0 , 0 ...} 


which is linear and isometric. Here J-+ denotes the bosonic Fock space over the orthogonal comple¬ 
ment of span{99}. We can think of U^ipN £ J^+ as describing the fluctuations around the condensate. 
We assume that the initial data 0jv £ exhibits complete condensation in 99 E L^(M^) in the 

sense that the fluctuations 0(0) ;= UipipN £ T+ can be bounded uniformly in N in the sense that for 
some C > 0 we have 

(0(o),dr(i-A)0(o)) <c 

®The symbol denotes the symmetrized tensor product, i. e. i/Ji (g)® • • • (g)s pN = (N!)“^ JDaeSiv ‘^'^( 1 ) ® ' ’ '®V’rr(N)- 
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independent of N. We let ipN evolve with the Hamiltonian 

N ^ N 

Hn = ~^xj + ^ V{xi — Xj). 
j=l i<j 

We know that 'i/jTv,* = exhibits condensation in the one-particle state (pt evolved according 

to the Hartree dynamics. To study the evolution of the fluctuations around the Hartree evolution we 
apply the map U^p^ : —)• defined analogously to on (Notice that the image 

space of depends on time because of the requirement of orthogonality to ipf) We find that 

||C/^,e-*^^V 7 V-</>(t)|| ^0 (iV^oo). (4.10) 

Here (j){t) G is the solution of the evolution equation 

idt(t){t) = Coo{t)cj){t) (4.11) 

with initial data ())(0) and a quadratic time-dependent generator Coo{t) very similar to (|4.ip (although 
not exactly the same because of the requirement of orthogonality). Eq. (14.101) shows that, by taking 
into account the quadratic dynamics of the fluctuations, we obtain a norm approximation for the 
full evolution In other words, writing ^(t) = (t),..., (t), 0, 0,... }, we find 

'lpN,t ^ 0s H-h 

with an error whose norm tends to zero as —)• oo. (It is also possible to check that the error is 

actually of order 

5 The Gross-Pitaevskii regime 

Another limit in which it is possible to approximate the many-body dynamics of bosonic systems by 
an effective one-particle equation is the Gross-Pitaevskii regime, which is relevant for the description 
of trapped Bose-Einstein condensates. At the microscopic level, a trapped Bose-Einstein condensate 
can be described as a gas of N bosons with Hamilton operator of the form 

N 

+ '^ext(^,)) +Y, N^ViNixi - Xj)). (5.1) 

j=l i<j 

Here lext is an external potential modelling the trap and the interaction is described by a smooth 
repulsive potential V (we have to assume that E > 0 pointwise) with short range (for convenience 
we will assume that V has compact support, although this is not really necessary). 

Scattering length. In (|5.1I) the interaction potential scales with the number of particles N so that 
its scattering length is of order N~^. Let us recall that the scattering length of a potential V is 
defined through the solution of the zero-energy scattering equation 

(-A + iE)/ = 0 (5.2) 

with the boundary condition f{x) —>• 1 for |x| —>■ oo. Under the assumption of compact support for 
V one can show that, for |a:| sufficiently large, 

fix) = 1 - —r 
\x\ 
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for an appropriate constant oq > 0, which is called the scattering length olV. It is a simple exercise 
to show that equivalently the scattering length oq can also be defined through the integral 


Sttuo 


V{x)f{x)dx 


(5.3) 


where again / denotes the solution of (|5.2p . From the point of view of physics, the scattering 
length Oq measures the effective range of the interaction potential; two quantum mechanical particles 
interacting through the potential V, when they are far apart, feel the other particle as a hard sphere 
with radius oq (in particular, the scattering length of a hard sphere potential coincides with the 
radius of the sphere). 

Notice that if oq denotes the scattering length of V, the scattering length of the rescaled potential 
N‘^V{N-) is given by a = a^/N. This follows because by simple scaling, from (|5.2p we find 

(-A + ^y(iv.)) /(iv-) = o 


and 


f{Nx) = 1 - 


Qq 

N\x 


= 1 - 


(ao/N) 

Ixl 


Ground state properties of trapped condensates. It was proven in [50] that the ground state 
energy per particle for the Hamiltonian (15.ip converges, as N ^ oo, towards the minimum of the 
Gross-Pitaevskii energy functional 

£gp{^) = j + Fextlv^l^ + 47rao|¥?|^] dx (5-4) 

over all one-particle wave functions (p G L^(M^) with ||(/?|| = 1. From (15.4p we conclude that in first 
approximation the ground state energy of the boson gas depends only on the scattering length of the 
interaction potential, not on its precise profile. 

It was then shown in [JH] that the ground state of exhibits complete condensation in the 

minimizer of the Gross-Pitaevskii energy functional (|5.4p . More precisely, the one-particle reduced 
density 7 ^^ associated with the ground state of was proven to satisfy 

^ l¥^)(</^l (N^oo). 

The interpretation of this result is straightforward: in the ground state of all particles, up to 

a fraction vanishing in the limit A —>■ 00 , are condensated in the one-particle state described by the 
unique minimizer of the Gross-Pitaevskii functional. 

Dynamics of initially trapped condensates. Since Gross-Pitaevskii theory has proved so successful 
in the description of the ground state properties of the Hamiltonian (|5.1I) . can it also be used to predict 
the time-evolution of initially trapped condensates? As in Section [ 2 ] we want to study the reaction of 
the system to a change of the external helds. At time t = 0 , we assume the boson gas to be prepared 
in the ground state ?/>Ar of the trapping Hamiltonian (15.ip . Then, the traps are switched off and the 
condensate starts to evolve by the translation invariant Hamiltonian 

N N 

= +^AV(iV(xi-x,)). (5.5) 

j=l i<j 

The next theorem, taken from [2711281 [29] . describes the resulting dynamics. (More recently, a similar 
statement has been shown in [54j.f 
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Theorem 5.1. Let V > 0 be a spherically symmetric, short range, bounded potential with scattering 
length qq. Consider a sequence ipN G satisfying 

• Finite energy per particle: there exists C > 0 such that {ipN, < CN; 

• Condensation: the one-particle reduced density 7 ^^ associated with 'ipN is such that 

^7^^^ ^ l<^)(</7| (N^oo) 

for a ip £ 

Let and let 7 ^ ^ be the one-particle reduced density associated with 'ipN,t- Then, for 

every fixed t G M, 

\Tt){v’t\ (5.6) 

as N ^ 00 , where ipt is the solution of the time-dependent Gross-Pitaevskii equation 

idtpt = -^Pt + ^T^ao\pt\^pt (5.7) 


with the initial condition po = p. 

Remarks. Convergence in (j5.6|) holds for example in the trace norm. Moreover (|5.6I) also im- 

(k) 

plies convergence of higher order reduced densities. If 7 )^) denotes the /c-particle reduced density 
associated with it is shown in [251 [ 22 ] that 

^ \Pt){Pt\'^’' 

as N ^ 00, for any hxed /c G N and t G K. 

Comparison with mean-field regime. Let us discuss the relation between the Gross-Pitaevskii 
limit, characterized by an interaction potential with scattering length of the order N~^, and the 
mean-field regime discussed in Section [21 Writing the interaction potential in (15.5p as 

N'^V{Nx) = ^vn{x) 

with vj\f{x) = N^V{Nx), one could think of the Gross-Pitaevskii limit as a mean-field limit with 
a potential vn that converges towards a delta distribution in the limit of large N. However, this 
interpretation of (|5.5p is quite misleading. Formally we have 

vn{x) = N^V{Nx) —)• bod{x) 

with bo = fV(x)dx. Hence, based on the mean-held interpretation of (|5.5p . we should expect the 
many-body evolution to be approximated, in the limit of large N, by the nonlinear Schrodinger 
equation 

idtpt = -Apt + bo{d * \pt\'^)pt = -Apt + bo\pt\‘^pt. 

This equation has the same form as the Gross-Pitaevskii equation (|5.7p but a different constant in 
front of the nonlinearity. The reason why the mean-held interpretation leads to a wrong coupling 
constant is that physically the two regimes are very different. While the mean-held regime is char¬ 
acterized by a large number of very weak collisions, in the Gross-Pitaevskii limit particles interact 
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very rarely (only when they are at distances of the order which is much smaller than the typ¬ 
ical distance among the particles) but when they do interact, the collisions are very strong. 

Because of these rare and strong collisions, the solution of the Schrodinger equation generated 
by the Hamiltonian (|5.5p develops a singular correlation structure, varying on the length scale 
which is then responsible for the emergence of the scattering length in the Gross-Pitaevskii equation 
(1521). Correlations among the particles therefore play a crucial role in the Gross-Pitaevskii regime, 
while they are negligible in the mean-field limit. 

Correlation structure. Let us now try to explain how the correlation structure affects the dynamics 
of the condensate. (Recently the correlation structure developed by the solution of the Schrodinger 
equation in the Gross-Pitaevskii limit has been studied in [20].) Let us normalize the one- and 
two-particle reduced densities associated with to ^ and y^^- They 

satisfy the differential equation 


+ (^ - l)Tr2 N'^ViNixi - X2)),y]vl 


;(i) 


( 2 ) 


(5.8) 


which is the first of the N coupled equations forming the BBGKY hierarchy (similar to ()2.9h in 
the mean-field setting). By assumption, at time t = 0, “ipNfi exhibits condensation, meaning that 
y^j^o for a G L^(R^). If condensation is approximately preserved by the time-evolution, 

we should expect that also for t 7 ^ 0 

(5.9) 

As for the two-particle reduced density, again we expect approximate factorization. Here, however, 
we should also take into account correlations. Assuming that correlations can be described by means 
of the solution of the zero-energy scattering equation (|5.2I) . we expect that the integral kernel of y)^") 
can be approximated by 

77 v,t(^i’^ 2 ;yi,y 2 ) fiN{xi -X 2 ))f{N{yi - y 2 ))(pt{xi)y:>t{x 2 )ipt{yi)^t{y 2 )- (5.10) 

Using this ansatz for y)^'), we find that the second term on the r. h. s. of (|5.8|) approximately has the 
integral kernel 


{N - 1 )TV 2 [YV(iV(xi - X2)),yg] {x;y) 

= {N-l) I dx2 {N^V{N{x - X2)) - N^V{N{y - X2))) yS^J(x, X2; y, X2) 

~ J dx2 {N^V{N{x - X2)) - N^V{N{y - X2))) ipt{x)Tp^{y)\ipt{x2)\'^ 

X f{N{x - X 2 ))f{N{y - X 2 )) 

~ 87rao (|¥’t(x)P - |</?t(y)|^) ipt{x)Tpt{y) 


(5.11) 


where we used ()5.3p and the fact that f{Nx) —?■ 1 weakly as N —>■ 00 . Inserting also (|5.9I] in the 
term on the 1. h. s. as well as in the first term on the r. h. s. of (15.8|) . we see that the Gross-Pitaevskii 
equation (|5.7p arises exactly as the self-consistent equation for cpt. Notice that the presence of the 
solution of the zero-energy scattering equation f{N-) in the ansatz ()5.10p for the two particle reduced 
density does not contradict the fact that, as N ^ 00 , y^^^ the correlation structure 

in (15.mil is non-trivial only for \xi — X 2 I < N~^ or |yi — ^ 2 ! ^ and it disappears in the limit 

N ^ 00 . Nevertheless it plays a crucial role in (15.1111 because it is multiplied with a very singular 
potential, varying on the same short scale. 
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Energy estimate. Since the presence of the solution f{N-) of the zero energy scattering equation 
in (IS.lOp plays such an important role, a rigorous derivation of the Gross-Pitaevskii equation (15.71) 
requires a proof that the solution of the many-body Schrodinger equation really develops a 
correlation structure and that, in good approximation, this correlation structure can be described 
by To reach this goal, it is useful to prove certain energy estimates, bounding appropriate 

Sobolev norms of tpN,t by moments of the Hamiltonian. 

To give an example of the analysis involved in this step, we prove a simple energy estimate for the 
case of small interaction potentials. More precisely, we define the following dimensionless quantity 
to measure the strength of the radial interaction V: 

POO 

yO = supr^l/(r) + / drrV{r) 
r>0 JQ 


where we introduced the short hand notation r = |x| and V{r) = V{x) if |x| = r. Smallness of p also 
implies that the solution / of the zero-energy scattering equation (15.2|) remains close to 1 (recall the 
boundary condition /(x) —?• 1 as |x| —>■ oo). In fact, one can show [271 Lemma D.l] that there exists 
a constant c > 0 with 

l-cp<f{x)<l, |v/(x)|<c^, |v2/(x)|<c^ (5.12) 

|x| |xp 

for all X G M^. Using (|5.12p . we can prove the following proposition, which is taken from [28j . 

Proposition 5.2. Let V > 0 be a spherically symmetric, short range, bounded potential with p > 0 
small enough. Let f denote the solution of the zero-energy scattering equation 115. S\) . Then 


> — dx 


V V __ 

"V(iV(xi-X2)) 


(5.13) 


for every ipN £ L‘^{'R^^). Here we used the eonvention x = (xi,... ,xn) G . 

Remark. Let be the solution of the Schrodinger equation with initial data 

Assume additionally (with respect to the assumptions of Theorem 15.ip that (V’at, < CiV^; 

this can be achieved with an approximation argument that we skip. Then (I5.13P implies that 


Ar2 



^ X2 


f{N{xi - X 2 )) 


2 

< < 2CN‘^ 


and therefore that 





'0Af,t(x) 

f{N{xi - X 2 )) 


2 

< C 


(5.14) 


uniformly in N and in t. Note that in (j5.14l) it is very important that we divide 4>N,t by /(A^(xi —X 2 )) 
before we take the two derivatives. Keeping in mind that /(x) = 1 — ao/|x| for |x| large enough, it 
is easy to check that 

J dx|v2(/(iVx))p ~ N 

and therefore we can expect that 

j dx\Vxi'^x2fjN,t{x)\'^ ^ N. 
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Only if we first divide V'TV.t by f{N{xi —X 2 )), removing in this way the singular correlation structure 
in the xi — X 2 variable, we obtain in (|5.14p a bound of order one. In other words, (I5.14p implies that, 
when \xi — X 2 I < the evolved wave function can be approximated by f{N{xi — X 2 )) times 
a function 'i/^jv,i(x)//(A^(xi — X 2 )) which varies on a larger length-scale. In this sense, we can say 
that the estimate p5.14p proves that has a short scale correlation structure, which can in good 
approximation be described by the function f(N(xi — X 2 ))- For this reason the bound (|5.14l) plays 
a crucial role in the proof of Theorem 15.11 (for the case of weak potential treated in [28]; for large 
potentials, a different energy estimate has been proved and applied in [29]i. 

Proof of Proposition \5.^ We write 

^ 1 
HN = '^hj with hj =+-'^N‘^V{N{xi - Xj)). 

i=i i¥=j 

Then, using the permutation symmetry of we hnd 

= N{N - l){ipN,hih2ipN) + N{il;N,hlijjN) > N{N - l){ipN,hih2ipN)- 

The positivity of the potential and the fact that for all j >3 

N^V{N{x, - x,))(-A,J = Vl^N^V{N{xi - Xj))V,, 


imply that 


N' 


{iPn, > N{N - -A^, + —V{N{xi - X 2 )) -A,, + —V{N{xi - X 2 )) M 




Next we define 4>n{^) = f{Af{xi — X 2 )) and write 

^ 2A(V/)(iV(xi-X2)) ^ , , N^{-Af){N{xi-X 2 )) 


f{N{xi - X 2 )) 
Hence 


— ^xi4^N 


fiN{xi - X2)) 


^ xi4>N + 


fiN{xi - X 2 )) 




-Aa;j + —V{N{xi - X 2 )) ) V’iV = fiN{xi - X 2 ))Li(j)]S[ 


(5.15) 

(5.16) 

(5.17) 


where we defined the operator 


Li 


-Axi4>N 


2N{Vf){N{xi-X2)) 

f{N{xi-X 2 )) 


The important observation is that the last term on the r. h. s. of p5.16p cancels exactly with the 
potential, because / satisfies the zero-energy scattering equation. Analogously, we find that 


+ ^V{N{xi - 2 : 2 ))^ 


V’iV = f{N{xi - X2))L2(I)N 


with 


L 2 


—Ax24>n + 


2iV(V/)(jV(xi-X2)) 

f{N{xi-X 2 )) 
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It is important to observe that Li,L 2 are the Laplace operator w. r.t. the weight In other 

words 


j - X2)(LiXAr)(x)^Ar(x) 


j dx/^(xi - X2)XAr(x)(-Li^Ar)(x) 


= j dxp{N{xi - X2))V2;iXAr(x)V3;i^Ar(x) 


and analogously for L 2 . 

From (I5.15P we conclude that 


{'iPn,HniPn) > N{N - 1) y dx/^(A^(xi - X2))Li(/)Ar(x)L2(/)Ar(x) 

= N{N - 1) y dy.f{N{xi - X 2 ))VxiL24>N{yd) 

= N{N-l) I d^f{N{xi-X2)) |V,,V,,</.Ar(x)|2 

+ iV(iV-I) y dx/2(Ar(xi-X2))V;;^(x)[V,,,L2]</.7v(x) 

= N{N-l) I d^f{N{xi-X2)) |V,2Vxi07v(x)|" 

+iv(iv-i) y dx/^(iv(xi - X2)) 

(5.18) 


In the first term on the r. h. s. we bound /^(iV(xi — X 2 )) > (1 — Cp)"^ from below (see (|5.12D i. As for 
the second term on the r. h. s. of (I5.18p . we notice that 

iV(V/)(iV(xi-X2)) ^ N^{V^f){N{xi-X2)) N^{Vff{Nixi-X2)) 

f{N{xi-X2)) f{N{xi-X2)) P{N{xi-X2)) 

The bounds (15.121) therefore imply that 


iV(V/)(iV(xi-X2)) 
f{N{xi - X 2 )) 


< Cp 

|Xl - X2 


for small p > 0. Hence 


iv(iv-i) y dx/^(iv(xi - X2)) ^x, 4 >n{^) 

< CpN{N - 1 ) [ (iX|- - -^|V3;i(?!)Ar(x)||Va;2(()Ar(x)| 

J \x 1 -x 2 r 

<CpN{N-l) y dx|V,,V,,(/.7v(x)|2 

and thus 

(V^TV, > (1 - Cp)N{N - 1) y dx \Vx,,VxAn{^)? . 

For p > 0 sufficiently small, we obtain the desired bound. □ 
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Strategy of the proof of Theorem \5.1[ As mentioned above the energy estimate (I5.13P and its 
corollary (15.141) play a crucial role in the proof of Theorem 15.11 because they can be used to identify 
the short-scale correlation structure characterizing the solution of the many-body Schrodinger 
equation. Let be the family of normalized reduced densities associated with satisfying 

= 1 for all iV, /c G N and t G M. Denote by a limit point of the sequence 

for N —)• oo. The estimate (j5.14p implies that, for large but hxed G N (along the appropriate 
subsequence), we can approximate 


k 

7 ^](xi,..., Xfc; yi,..., y/c) ^ f{N{xi - xj)) 7 ^|i(xi,..., yi,..., y/c) (5.19) 

i<j 

This allows us to derive, starting from the BBGKY hierarchy for the reduced densities 
the inhnite hierarchy 


(fc) _ 




K 

n 

i=i 


-A, 


~{k) 


+ Svroo ^ Trfc+i 


6{xj - Xfc+i),7, 


(k)- 

oo,t 


(5.20) 


for the limit point { 7 ^\}fc>i. The presence of the correlation structure in ()5.19l) is the reason why the 
scattering length oq arises in the infinite hierarchy (|5.20p . What is still missing to show Theorem 15.11 
is a proof of the uniqueness of the solution of the infinite hierarchy (|5.2np . Before proving uniqueness, 
it is important to understand in which space of families of densities { 7 ^|i}fc>i uniqueness should be 
shown; proving uniqueness in a smaller space is of course easier, but it also requires to show that 
every limit point of the sequence is contained in that small space. It turns out that a 

possible choice of the space of densities is the Sobolev-type space 

ni = {{7^''^}fc>i: 30 0 Vfe G N Tr (1 - ... (1 - A^JyW < C'^}. 

Hence, to conclude the proof of Theorem l5.ll one has to prove that every limit point of { 7 ^j}a:>i is 
in the space Tii, and further one has to show that, for any initial condition 7 ^^o ~ ^ 

there is at most one solution { 7 ^\}fc>i of the inhnite hierarchy (|5.2UD which is in Tii for all t G M. 

The proof that every limit point {^^\}k>i of the sequence satishes the a-priori bounds 

Tr(l-A,j...(l-A,j7gi (5-21) 

for all /c > 1 is based on energy estimates similar to Proposition 15.21 but involving also higher 
moments of Hj\f. A challenge is the fact that for hnite N, the reduced densities 7 )^') cannot satisfy 
(I5.2ip . at least not uniformly in N, because of the presence of the correlation structure (taking 
derivatives of the correlation functions f{N(xi — Xj)) produces factors of N). Only after taking the 
limit A — 00 , the correlation structure disappears and one can expect (I5.2ip to hold true. To obtain 
(|5.2ip . one needs therefore to show estimates of the form 


j dxeeix) |Va;i. 


..V,,V’iv,i(x)P <0^^ 


(5.22) 


where the cutoff &£ satishes &i(x) ~ 0 if there exists i < k and j < A with \xi — Xj\ < i, while 
©^(x) ~ 1 if |xj — Xj\ ^ i for all z < A; and j G {1,.. ., A} with j / z. If £ is sufficiently large (it 
turns out that one needs A£^ ^ 1), the cutoff removes all singularities of due to 
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the correlation structure (because correlations are only important when particles are close to each 
other). At the same time, if N(.^ <C 1, the effect of the cutoff turns out to be negligible in the limit 
of large N; hence ()5.21l) follows from ()5.22D . choosing <C ^ <C . 

The proof of the uniqueness of the infinite hierarchy (|5.2Up in the class T-Li, given in [2Tj, is 
based on an diagrammatic expansion in Feynman graphs. Since the singularity of the interaction 
(the 5-function (I5.20p l cannot be controlled by the kinetic energy (the "Hi-norm), one needs to make 
use of the dispersive properties of the free evolution generated by the hrst term on the r. h. s. of 
(I5.20p . To perform the analysis it is convenient to switch to Fourier space. Thanks to the decay in 
momentum characterizing every family of densities vnT-Li and to the decay of the propagators of the 
free evolution, the contribution associated to every Feynman graph is convergent in the ultraviolet 
regime. Further details can be found in m- 

A different and shorter approach to prove the uniqueness of the infinite hierarchy ()2.10p was 
later proposed in [33] based on certain space-time estimates for the densities Recently [22] . 

this approach was applied to deduce uniqueness in the space "Hi; an important ingredient was the 
quantum de Finetti theorem. 

Coherent states approach in the Gross-Pitaevskii regime. Theorem 15.11 proves the convergence 
towards the Gross-Pitaevskii dynamics without control on the rate of the convergence. Since in 
real systems the number of particles N is large but finite (e. g. ~ 1000 in very dilute samples of 
Bose-Einstein condensates), it is important to know how large N must be in order for the Gross- 
Pitaevskii equation to become a good approximation of the many-body quantum dynamics. Can 
the coherent states approach presented in Section [3] for the mean-field regime also be applied in the 
Gross-Pitaevskii limit to obtain an explicit bound on the error? 

Let us hrst try to proceed naively, following exactly the same strategy as in the mean-held case. 
We switch to the bosonic Fock space J- = 0^>qL^(M^”) on which we dehne the Hamilton operator 

T-Ln = J dxVxalVxax + ^ J dxdyN‘^V{N{x - y))ala*yayax. (5.23) 

We consider a coherent state initial data with a G L^(R^) with ||(/?|| = 1. The 

expected number of particles in the state 'kjv is N. We let 'kjv evolve and try to approximate it with 
a new coherent state. To this end we dehne the huctuation vector G by 

e-*^^'W(ViV(/p)L! = W{VNipt)iN,u 


which can be rewritten as 

iN,t = WAr(t;0)n 

with the huctuation dynamics 

U^it- s) = W{y/Nipt)*e-^'^^^W{y/N^s)- (5.24) 

Instead of choosing ipt to solve the Gross-Pitaevskii equation (15.7p . it is convenient to consider the 
solution of the modihed Gross-Pitaevskii equatiorj^ 

idm = + {N^V{N-)f{N-) * (5.25) 

®Of course, this equation depends on N, and so do its solutions. So when we talk about convergence of this equation 
to the Gross-Pitaevskii equation we mean convergence of its A-dependent solutions to solutions of the Gross-Pitaevskii 
equation. Despite solutions being A-dependent we don’t put an extra A-index to keep the notation light. 
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with the same initial data tpQ = tp. As —)• oo, 

N^V{Nx)f{Nx) —)• 87rao6{x) 

and the solntion of the modified eqnation (15.251) can be easily shown to converge towards the solntion 
of the Gross-Pitaevskii equation (15.7p . with the rate N~^ (at least for sufficiently regular initial data). 

Similarly as in the mean-field case, to show that {^/Nip)Q can be approximated by 

the evolved coherent state it is enough to prove that the fluctuation vector stays close to the 
vacuum. More precisely, it is enough to control the growth of the expectation of the number of 
particles operator 

{l(N{t;0)n,MUNit;0)n) (5.26) 

with respect to the fluctuation dynamics To this end, we compute the (time-dependent) 

generator C^it) of Z^Ar(t;s), defined by the equation 

idtUNit; s) = CN{t)VlN{t; s), Un{s; s) = 1 


and thus given by 


CN{t) 


idtW*{yfNLpt)] W{yfNpt) + W*{y/Npt)nNW{yfNpt). 


On the one hand we have 


idtW*{VNpt) 


W{VNipt) = Ci{N,t) - Vn [a*{idt(pt) + a{idtpt)] 


(5.27) 


for an unimportant constant Ci{N,t). On the other hand, using (|3.17p . we find 

w*{VNpt)nNW{VNipt) 

= C2{N,t) + ^/iV [a*{-A^ipt + (iW(iV-) * \ipt\^)pt) + a{-A,ipt + {N^V{N-) * \pt\^)ipt)] 

+ j dxVxalVxax + J dx{N^V{N-) * \ipt\‘^){x)alax 

+ j dxdyN^V{N{x-y))pt{x)Tp^{y)alay 

+ ^j dxdy N^V{N{x - y)) [pt{x)n{y)alal+ Tp^{x)Tp^{y)a^ay\ 

+ ^/ dxdyN^V{N{x-y))ala*yayax 

(5.28) 


where C 2 {N,t) denotes another unimportant constant. Recall that in the mean-field regime the 
linear terms (linear in creation and annihilation operators) arising from the two contributions to the 
generator (|3.26p canceled exactly due to the choice of pt as a solution of the Hartree equation. Here 
in the Gross-Pitaevskii regime, the cancellation is not complete] in fact the creation and annihilation 
operators on the r. h. s. of (|5.27p have the argument 

idm = -Apt + {N^V{N-)f{N-) * \ptf)pt 
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while the argument of the creation and annihilation operators appearing in the linear terms on the 
r. h. s. of (|5.28l) does not contain the solution f{N-) of the zero-energy scattering equation (I5.2|) . It 
follows that 


C^it) = C{N,t) + Vn [a*{{N^V{N-)u;{N-) * + a{{N^V{N-)u;{N-) * 

+ J dxV ri;alV + j rf® ( F ( iV' ) * | ^ ) (x ) Oj; 

-f j dxdyN^V{N{x - y))^pt{x)Tp^{y)alay 

+ ^J dxdy N^V{N{x - y)) [ft{x)ipt{y)ala*y+ Tp^{x)Tp^{y)a^ay\ 

+ ^/ dxdyN^V{N{x - y))ala*yayax 


where we set a; = 1 — / (and where C{N,t) = Ci{N,t) + C 2 {N,t) is another unimportant constant). 
Hence, in the Gross-Pitaevskii regime, the generator of the fluctuation dynamics (15.2411 contains a 
contribution which is linear in creation and annihilation operators and therefore does not commute 
with the number of particles operator. This contribution is of order y/N (at least formally), and 
for this reason, it is impossible to bound the growth of the expectation of the number of particles 
operator (|5.26p uniformly in N. 

Modified fluctuation dynamics. The reason why this naive approach does not work is the fact 
that we are trying to approximate the many-body evolution of a coherent state initial data by 
an evolved coherent state, completely neglecting the correlation among the particles. (Recall that 
coherent states in each re-particle component are completely factorized and thus do not have any 
correlations.) Since we know correlations to be very important in the Gross-Pitaevskii regime, it 
should not surprise us that fluctuations with respect to the evolved coherent state are too large. 

To take into account the correlation structure developed by the time-evolution we define the 
kernel 

kfix] y) = -Nu}{N{x - y))(pt{x)ipt{y) (5.30) 

where, as before, ipt denotes the solution of the modified Gross-Pitaevskii equation (15.2511 and io = 
1 — / (where / is the solution of the zero-energy scattering equation (15.21) 1. Recall that lo{x) ~ ao/|x| 
for |x| S> 1, while in accordance with (j5.12p it is regular for |x| <C 1; it is useful to think of the function 
Nuj{N{x — y)) as ao/|x — y\ but regularized for \x — y\ < 1/N. Using the kernel kt we define the 
unitary operator Tt, acting on the Pock space by 

Tt = exp Q j dxdy (/ct(x; 7/)a*a* -ktix;y)axay)^ . (5.31) 

Since the exponent is quadratic in creation and annihilation operators, Tt implements a Bogoliubov 
transformation. It is even possible to explicitly compute its action on the creation and annihilation 
operators, yielding 


Tta{f)Tt = o(coshfc,/) -F a*(sinhfct/) 
Tta*{f)Tt = a*{coshkJ) + a(sinhfc,7) 


35 













for any / G Here we use the notation cosh^j and sinh^^ for the linear operators on 

given by 

coshfc, = = Z] + 

n>0 ^ n>0 ^ r 

where products of kt and kt have to be understood as products of operators (we identify the function 
kt G L^(M^ X M^) with the operator having kt as its integral kernel). We use the unitary operator Tt 
to approximate the correlation structure developed by the many-body evolution. We consider the 
evolution of initial data having the form 

'kTv = kh(v^</^)ro^7V (5.33) 

for a with only few particles (think of = O for simplicity), and we try to approximate its 
evolution with an evolved state of the same form. As discussed in Appendix!^ ()5.33p is a natural 
class of initial data, approximating the ground state of the Bose-Einstein condensate trapped in a 
volume of order one (the point here is that the Bogoliubov transformation Tq generates the correct 
correlation structure, which is crucial on the one hand to reach the ground state energy and, on the 
other hand, to follow the many-body dynamics). 

To approximate the evolution of (|5.33l) . we define the new fluctuation vector requiring that 

e-'^^'W(ViV(^)roe7v = W{y/N^t)TtiN,t. 

Equivalently 0)^ with the modified fluctuation dynamics 

UNit] s) = T:W*{VNipt)e-^^^^^-^^W{VNips)T,. (5.34) 

Let us compute the kernel of the one-particle reduced density 7 ^^ associated with = 


= {^N,t,a*yax'i>N,t) 

= {W{y/Nipt)Tt^N,t, a;a^W{VNipt)Tt^N,t) 

= {^N,t,Tt* [a; + VNlf^iy)) (a, + VNipt{x)) Tt^N,t) 


Hence, for any one-particle observable J on L^(]R^), we hnd 

Tr J - N\y^t){^t\) = VN{(N,t,T; [a*{J^t) + a{Jipt)] Tt^Ny) + {iN,t,T*dr{J)TtiN,t) 
and therefore 


Tr 


7 $;) - N\^t){^t\ < VN{^N,t, T*M^/^TtiN,t) + {iN,u Tt*AfmN,t). 


Using (I5.32p . we conclude that 
and therefore that 


Tr 


T^MTt < C{M + 1) 

7S;j - N\^t){‘Pt\ < CVN{CN,t, (AT + l)CiV7)- 


(5.35) 


(5.36) 


is interesting to note that, while the introduction of Tt changes the energy by a contribution of order N (see 
Appendix Eq. (16.2811 shows that the change in the number of particles is only of order one. 


36 
















This means that to get a bound on the rate of the convergence of the many-body evolution towards 
the Gross-Pitaevskii dynamics (proportional to ^/N), it is enough to control the growth of the number 
of particles with respect to the modified fluctuation dynamics, i. e. to control the expectation value 

(5.37) 

uniformly in N. 

Generator of modified fluctuation dynamics. To estimate the expectation (|5.37|1 we compute the 
generator of V(N{t] s), defined by the equation 

idtV(]\f{t] s) = CN{t)UN{t] s) with Un{s] s) = 1 


and therefore given by 

CN{t) = {idtTn Tt + (5.38) 

where Cisf{t) is the generator computed in (I5.29p . The first term on the r.h. s. of (|5.38|1 is an 
expression quadratic in creation and annihilation operators and can be bounded as 

±[{idtT;)Tt,Af]<C{Af + l). 


Hence, its contribution to the growth of the number of particles operator can be controlled similarly 
as we bounded (j3.29p in the mean-field case. Let us now focus on the second term on the r. h. s. 
of (j5.38D . As discussed between Eq.s (j5.28ll and (I5.29p . the generator C]\f{t) contains a large term 
(proportional to \//V) linear in creation and annihilation operators. After conjugation with Tt, this 
term is given by 

Vn j dxdyN^V{N{x - y)MN{x - y))\^t{y)\^Mx)T:alTt + h.c. (5.39) 

At the same time, the generator Cisf{t) (I5.29p contains a term which is cubic in creation and annihi¬ 
lation operators. After conjugating with Tt we find this contribution to be, introducing the notation 
ch.x{z) = coshfcj(z;x) and sha,(z) = sinh( 2 :;x), 


^/N . 


dxdyN^V{N{x - y))ipt{y)Tf alaya^Tt + h.c. 


1 

y/N 

y/N 


J dxdyN^V{N{x - y))Lpt{y)Tf alTt{a{c\iy) + a*(shy))(a(cha,) a*{shx)) -h h.c. 
dxdyN^V{N{x - y))ipt{y)Tf a^Tta*{sh.y){a{ch^) + a*(sha;)) 

1= J dxdyN^V{N{x-y))ipt{y)TfalTta{chy)a{ch^) 


(5.40) 


y/N 

1 


y/N . 

+ h.c. 


dxdyN^V{N{x - y))ipt{y)Tf alTta{diy)a*{shfij 


In the last summand, the operators a(chj^) and a*(sha,) are not in normal order (a product of creation 
and annihilation operators is said to be normal ordered if all creation operators are to the left of all 
annihilation operators). Putting them in normal order produces a commutator term which is linear 
in creation and annihilation operators. Since cosh^^ = 1 -|- 0{kt/ and sinh^j = kt + 0{kt), we find 


(cha;,shy) = j dzcoshkfiz-,x)smhkfiz;y) j dz6{x - z)kt{z-,y) = kt{x-,y) (5.41) 
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up to terms which are regular in the variable x — y (more precisely, higher powers of have kernels 
that are regular on the diagonal and therefore their contribution to (I5.4U|) can be shown to be 
negligible; in contrast, kt{x\y) behaves like ao|x —for \x — y\ ^ 1/-^)- With this approximation, 
we find from (|5.4Up 


J dxdyN^V{N{x - y))Tp^{y)T*alayaxTt+ h.c. 

~ I +h.c. 

~ -VlV [ dxdyN^V{N{x - y))uj{N{x - y))\ipt{y)\'^ipt{x)T^ a^Tt+£+ li.c. 


(5.42) 


where £ includes all terms proportional to T^a%Tt multiplied with a normal ordered quadratic ex¬ 
pression in a and a*. It is easy to check that all terms in £ after expanding T^a*Tt can be written in 
normal order up to negligible errors, and the normal ordered expression can be bounded using the 
number operator Af. In contrast, the main term on the r. h. s. of (15.421) can not be bounded; however 
it cancels exactly the large linear contribution in 15. S9\) . 

There is another important cancellation produced by the conjugation with the Bogolubov trans¬ 
formation Tt- On the one hand, we have the two quadratic contributions 


and 


J dxVxia*{ch^) + o(sh3;)Va,(a(cha,) -h a*(sha;)) 

= j dxdyVxal(yxkt){x,y)al + £i 

= J dxdy{Auj){N{x - y))ipt{x)ipt{y)al + £2 

dxdyN^V{N{x - y)) tpt{x)^pt{y) {a* {ch.^) + a(shj,))(a*(chy) -h a(shy)) 

= J dxdyN^V{N{x - y))ipt{x)ipt{y)alal + £3 


(5.43) 


(5.44) 


where £i ,£2 and £3 denote error terms, which can be controlled by the expectation of the number 
of particles operator Af and of the kinetic energy operator 


/■ 


/C = / dxVxofx^xO-x- 


1 


On the other hand, from the quartic term in (15.291) we find after conjugation with Tt 

J dxdy N^V{N{x — y)){a*{chx) + a(sh 3 ;)(a*(chj^) -|- a(shy))(a(chy) -|- a*(shy))(a(ch 2 ;) -|- a*(sh 3 ;)) 

= ^ J dxdyN^V{N{x - y))a*(ch 3 ,)a*(chy)a(chy)a*(sha;) + £^ 

= J dxdyN^V{N{x - y))kt{x-,y)a*x.al + £3 
= - j dxdyN^V{N{x - y))uj{N{x - y))(^t(x)<y9t(y)o>* -k £5 

(5.45) 
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where the error term 8 ^ can be bounded by the expectation of M, of M'^/N and of the quartic 
potential term 

J dxdyN^V{N{x — y))a*a*ayax ■ 

Combining the first terms of the r. h. s. of (Ih.d.fj) . (|5.44D and (|5.45D . we find 


/ dxdy 


-A + -V]{l-uj) 


{N{x - y))ift{x)<ft{y)ala* = 0 


(5.46) 


because we chose / = 1 — w as the solution of the zero-energy scattering equation (|5.2p . 

Taking into account these two important cancellations, the generator of the modified fluctuation 
dynamics (|5.34p can be controlled by the number of particles operator M, by /N and by the 
Hamiltonian (|5.23p . i. e. 

± CN{t) < (5.47) 

for some C, c > 0 (independent of N and t). The time dependence on the r. h. s. of the last equation 
arises through high Sobolev norms of the solution (pt of the modified Gross-Pitaevskii equation (15.251) . 

Growth of fluctuations in the Gross-Pitaevskii regime. The final step is to control the growth of 
the number of particles operator. In a very similar way as used to prove (15.471) . we also obtain the 
two bounds 


CN{t),N 


< Ce" 




[M+^ + CNit)^ 


(5.48) 


Furthermore we have a simple bound where the number of fluctuations is just bounded by the total 
number of particles, the ‘worst case’: 

(7/Ar(t;0)^Ar, —7/Ar(t;0)^Ar) < {UN{t]t))^N ( Cat, —Ctv) 

With these ingredients the next proposition can be proven by Gronwall’s Lemma. 

Proposition 5.3. Let H > 0, spherically symmetric and short range. Let G T, with 

+ {f,N,PNf,N) < C 


for some C > 0 (independent of N). Then there exists constants K,c> 0 (independent of N and t) 
such that 


(7/(t;0)CAf,AA7/(t;0)^Ar) < iFexp(exp(c|t|)) 


(5.49) 


for all t G K. 


The fast deterioration in time of the bound (|5.49l) is a consequence of the time-dependence in 
(I5.47P and (I5.48p . arising from bounds for high Sobolev norms of ipf If one assumes that 1199*11^4 < C 
uniformly in time, the bound (j5.49D becomes a simple exponential. 

Quantitative convergence towards the Gross-Pitaevskii equation. From (I5.36P and Proposition 
15.31 we obtain the following theorem, which has been proved in |12] . 
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Theorem 5.4. Let V > 0, spherically symmetric and short range. Let ip G and G 

such that 

+ {CN,hlNCN) < C 

for some constant C > 0 (independent of N). Consider the family of initial data 

= W{VNp)ToCn 


and denote by 7 ^^ the one-particle reduced density associated with the evolution = e 
0 / 4 ^TV- Then 


Tr 


7tvt - ^\Tt){v>t\ < exp(exp(c|t|)) 


for all t G M. 


6 Mean-field regime for fermionic systems 

So far we have discussed the time evolution of systems of interacting bosons. In this section, we will 
focus instead on fermionic systems. We will consider systems of N fermions initially confined to a 
volume of order one by a suitable external potential, e. g. an electromagnetic trap. We are interested 
in the evolution of such systems resulting from a change of the external field, e. g. switching off the 
trap. 

The Hamilton operator of the trapped N fermion gas has the form 

N N 

^7” = + ^ext(x,)) + A V{Xi - Xj) . (6.1) 

j=l i<j 

It acts on the subspace of consisting of functions which are antisymmetric w. r. t. 

permutation of the N particles. The interaction potential V varies on the same length scale as the 
one characterizing the confining potential text- For this reason, each particle interacts with all the 
remaining (N — 1) particles, producing a potential energy of the order AiV^. The mean-field regime 
is defined by choosing the coupling constant A so that the potential energy is typically of the same 
size as the kinetic energy. While the potential energy is independent of the statistics of the particles, 
the antisymmetry of the wave function plays an important role for the kinetic energy. In order to 
understand this point it is instructive to consider the following simple example. 

The free Fermi gas. Consider a gas of N non-interacting fermions in the three-dimensional torus 
T^. The Hamiltonian of the system is 




acting on 
determinants: 


2^2('j’3Af)_ eigenstates of can be computed explicitly; they are given by Slater 

nts: 


V^Slater(^1 ? • • • ; ^7v) v/iVT ^ ^ fl(,^n(^l)^ ' ' ' fN(,^n{N)^ ? 

TrES'iv 


where 

fi{x) = fpi{x) = e^P^'^, pi^2-Kl? . 
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For an arbitrary choice of N pairwise distinct momenta pi,... ,p]\f G 2tx1?^ the energy of the cor¬ 
responding Slater determinant is given by Thus, the ground state of is obtained by 

choosing N pairwise distinct vectors in 27rZ^ minimizing The condition that all 

momenta must be distinct follows from the observation that the Slater determinant vanishes if two 
momenta coincide. This is an expression of the Pauli exclusion principle, stating that there cannot 
be two fermions in the same one-particle state. It follows that the ground state energy of is 

approximated by filling the Fermi ball {p G 27 rZ^, \p\ < where c = 27 r/(^)^/^. We have (up 

to errors of lower order) 


inf spec(Ff^‘^^) ~ 


E 

\p\<cN^/^ 



for N ^ 1. 


Hence the ground state energy, which is purely kinetic in this simple example, is of the order 
and thus much larger than N. 

Fermionic mean-field regime. The fact that the kinetic energy of a Fermi gas scales as is a 
general property of N fermion systems trapped in a volume of order one. This is made rigorous by 
the Lieb-Thirring kinetic energy inequality, which states that for any ip G 


ifi’, 


N 

i=i 


p]l^{x)dx 


( 6 . 2 ) 


for a universal constant C. Here we defined the density p.^ associated with ip through 

p^{x) = N j dx2--- dxN IpJ^x, X2 ,..., xn)\‘^ ■ (6.3) 

Applying (j6.2p to the ground state ipN oi we conclude that, for well-behaved densities p.^, 


N 

{ipN, .iPn) > (6.4) 

i=l 

Therefore, in order to have a nontrivial theory in the limit N —>■ oo, we are led to choose the coupling 
constant A in Eq. m as X = N In the same spirit, the external potential should also be scaled 
with N: taking I4xt(3^) = N‘^^^Vext{x) we make sure that all components of the Hamiltonian scale as 
This completes the definition of the fermionic mean-field regime. Notice that if 

..Jo if X G A 

extfa^j ( oo otherwise 


for a domain A G the rescaling of the external potential is trivial. 

Hartree-Fock theory. In this regime, the ground state of the system is expected to be well 
approximated by a Slater determinant ipsiatevixi, • • • ,xn) = det{fi{xj)) with a suitable choice 

of the orthonormal system in T^(M^). The one-particle reduced density associated with V'siater 

is given by 

N 

OJN = A^Tr2,...,Ar|V'siater)(^/’Slater| = 

i=i 

Notice that con is the orthogonal projection onto the subspace of L^(M^) spanned by the orbitals 
{fj}f=i- Slater determinants are quasi-free states; that is, they are completely characterized by their 
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one-particle reduced density lon in the sense that all higher order densities can be computed through 
UN using Wick’s rule (|6.4UI) . In particular, the energy of a Slater determinant can be expressed in 
terms of only un, through the Hartree-Fock energy functional 


£ hf {’^ n ) = (V’Slater,-f^^^^V’Slater) 

= Tr(-A + I4xt)w + 

Hence, we expect the one-particle reduced density of the ground state of (16.ip to be well approxi¬ 
mated, for A = and N large, by the minimizer of (16.61) among all orthogonal projections un 

on L^(]R^) with TrcjAr = N. 

Thomas-Fermi theory. The Hartree-Fock energy functional (|6.6I) still depends on N. So what 
happens if we take the limit N ^ oo here? We obtain the next, coarser, degree of approximation 
which is Thomas-Fermi theory. More precisely, after dividing by the Hartree-Fock ground 

state energy is expected to be close to the minimum of the Thomas-Fermi energy functional 


j dxdyV{x - y)[u{x-,x)u{y,y) - |a;(x;?/)p 


( 6 . 6 ) 


Srpip) = ^CTF j dxp{xf/^ + 


dxVey±{x)p{x) + - 


dxdyV{x - y)p{x)p{y) 


(6.7) 


over all densities p E L®/^(R^) with p > 0 and ||/9||i = 1. The minimizer ptf satishes the Thomas- 
Fermi equation 


ctfPtf{x) 


{p - (l)TF{x)fl‘^ 


( 6 . 8 ) 


with the Thomas-Fermi potential 4>tf = Vext + V * ptf generated by the trap and by the self- 
consistent interaction produced by ptf- In (|6.8I) the parameter p plays the role of the chemical 
potential; it has be chosen so that ||/?tf||i = 1- 

Notice that ptf should be interpreted as the (normalized) configuration space density of the 
fermions in the ground state. With it, we can approximate the minimizer of the Hartree-Fock energy 
functional (|6.6I) by the Weyl quantization 

Wiv(x;y) = J dvM (6.9) 

of the characteristic function M(x,v) = x{\x\ < cPtf^^)) constant c > 0 is chosen so that 
Tiun = FI and is independent of A) 1^. Physically (|6.9n corresponds to the idea that in analogy to 
the case of free fermions, the ground state of (|6.6p (and of the corresponding many-body Hamiltonian 
(j6.ip l can be approximated by filling the Fermi ball; here, however, this procedure is implemented 
locally, with the number of particles in the local Fermi ball dictated by the Thomas-Fermi density 
Ptf- 

Dynamics. We consider initial data close to the trapped ground state of (16.111 and its evolution 
resulting from a change of the external potential. For simplicity we assume that at time t = 0 
the external traps are just switched off. The subsequent evolution is governed by the Schrodinger 
equation 

idr'lpN^ = HniPn,t , IpNfi = 'IpN 

with 

^ 1 

Hjv = ^-A„, + -Xj) . (6.10) 

i=l i<j 

^^However wn is only approximately an projection. 
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We now identify the relevant time scale on which the system undergoes macroscopic changes. From 
the discussion above, we expect the kinetic energy per particle to be of order hence, the classical 

velocity of the particles is typically of order Therefore the natural time scale of the evolution is 

of order After rescaling time by introducing the variable t = the Schrodinger equation 

takes the form 

N 

iN^/^dt'iljN,t= ■ (6.11) 

2=1 i<j 

Setting e = and multiplying the 1. h. s. and the r. h. s. of (16.111) by we find 

N ^ 

iedt'ilJN,t = ( ■ ( 6 . 12 ) 

2=1 i<j 

In (I6.12P we recover the N~^ coupling constant characterizing the mean-field limit of bosonic systems. 
Here, however, the mean-field limit is naturally linked with a semiclassical limit where e = 
plays the role of Planck’s constant. 

We are interested in the evolution (j6.12l) for initial data approximating Slater determinants with 
reduced density lon minimizing (|6.6|) . We expect that 'ipN,t stays close to a Slater determinant whose 
one-particle reduced density evolves according to the time-dependent Hartree-Fock equation 

iedtU}N,t = [hHF{t),ujN,t] , with hnFit) = -e^A + V * pt- Xt , (6.13) 


and with initial data = ^N- Here, [A, B] = AB — BA denotes the commutator of the operators 
A and B. Moreover, pt{x) = N~^LJN,t{x] x) is the configuration space density associated with 
the convolution V * pt, called the direct term, acts as a multiplication operator on L^(M^) and Xt is 
the exchange term, an operator with integral kernel given by Xt{x;y) = N~^V{x — y)u}]\[,t{x',y)- 
Like the Hartree-Fock energy functional, also the time-dependent Hartree-Fock equation (|6.13|) 
still depends on A^. As A" —oo, it is expected to converge towards the classical Vlasov equation. 
More precisely, let us define the Wigner transform 


WN,t{x,v) 



(6.14) 


of the reduced density As A —)• oo, is expected to converge (in an appropriate sense) to 

a probability density Woo^t on phase space which satishes the Vlasov equation 

{dt + v- V,,)Woo,i(x, v) = {V^V * pf){x) • V,Woo,i(x, v) , 

where p''^{x) = f dvWoo,t{x,v) is the configuration space density associated to Woo,t- The conver¬ 
gence of the Hartree-Fock (or Hartrej^ dynamics towards the Vlasov evolution has been proved to 
hold in a weak topology for sufficiently regular initial data H7v,o and for a large class of potentials 
including the Coulomb potential in m- More recently, precise bounds on the rate of the convergence 
from the Hartree equation to the Vlasov equation for regular initial data and for regular potentials 
have been established in mm- Notice, however, that in all these works the assumptions on the 
initial data exclude the case that Wn,o{x,v) ~ x(|x| < cp^/^(x)) for any reasonable density p, which 
would be interesting to consider for equilibrium states at zero temperature. 

^■^The exchange term Xt can be neglected with an error smaller than that of the step from many-body quantum 
mechanics to Hartree-Fock, see m- 
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Semiclassical structure. To establish the convergence of the many-body evolution to the time- 
dependent Hartree-Fock equation (|6.13l) . it is very important to take into account the structure of 
the initial Slater determinant. We are particularly interested in Slater determinants whose reduced 
density un is chosen close to the minimizer of the Hartree-Fock energy (16.61) . In this situation we 
expect the integral kernel ojm{x\ y) to decay to zero for |x — ?/| S> e. At the same time, we expect it 
to be regular and to vary on scales of order one in the x + y direction. 

To explain this point better, let us consider again the free Fermi gas on the torus. In this case, 
the one-particle reduced density of the ground state is given by 

0 JN{x\y) = g*P-(a:-y) ~ iV / dp = Nip(-—as N ^ oo, (6.15) 


for some function ip which can be explicitly calculated. As expected, cuw varies on the short scale e 
in the x — y variable. In this elementary case the (x -|- y)-dependence is trivial because Vext{x) = 0 
in the box A. 

The formula (16.9p suggests that the approximation wat of the minimizer of the Hartree-Fock energy 
functional dST]) for an interacting system in the mean-held regime exhibits the same semiclassical 
structure; it varies on the scale e in the (x — y)-variable and on a scale of order one in its (x + y)- 
dependence. 

A convenient way to quantify the separation of scales of a ;at is given by estimates on the commu¬ 
tators 

[a;Ar,a^] and ■ (6.16) 

The integral kernel of [ujntx] is given by 

> x]ix-, y) = (x- y)uN{x; y) = J dv (V^M) 



while 


[a;Ar,eV](x;y) = -e + '^y) u:N{x;y) 




x — y 
£ 


1 /3 

Hence, taking M{x,v) = xi\p\ < cpr^pix)), semiclassical analysis suggests that 


Tr |[ti;Ar,a:]| ^CNe J dxdv\{VyM){x,v)\ <CNe J dxppp{x) <CNe (6-17) 

and 

Tr |[wiv,eV]| < CNe J dx\VpTF{x)\ < CNe (6.18) 

under reasonable assumptions on V and 14xt (to guarantee that ptf & L^/^(M^) and Vptf ^ L^(M^) 
with iV-independent norms). 

In the following we will study the time evolution of approximate Slater determinants with one- 
particle reduced density satisfying the commutator estimates (I6.17p and (I6.18p . This assumption 
is motivated by the arguments given above, since we are interested in the evolution of initial data 
approximating the ground state of in the mean-held regime. 

Rigorous results. The mathematical literature concerning the derivation of effective evolution 
equations for interacting fermions is much smaller than for the bosonic case. In the mean-held limit 
we are interested in, the hrst derivation of the Vlasov equation from quantum many-body dynamics 
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has been given in [52] for analytic interaction potentials. Shortly after, this result has been extended 
in [62| to twice differentiable potentials. The first derivation of the Hartree-Fock equation has been 
given much later in [25|, for analytic interaction potentials and short times. 

In the rest of this section we will present the results of where the convergence towards the 
Hartree-Fock equation has been established with control of the rate of convergence for a much larger 
class of potentials and for all times (notice that the results of [T5| were also extended to the case of 
pseudorelativistic particles in (Ej). 

Concerning other scaling limits, the Hartree-Fock equation has been derived in |10] and (for a 
Coulomb potential) [3T] starting from a many-body Hamilton operator of the form (|6.ip . but with 
A = N~^ (in this case, there is no link to a semiclassical limit). Recently, convergence in presence of 
a regularized Coulomb interaction with a coupling constant has been given in |53j . 

Fock space. In order to discuss the results of [15] we need to introduce the fermionic Fock space 
over L^(M^), which is defined as the direct sum 

^ = C©0L2(M3n^ 
n>l 


where denotes the subspace of L^(M^"') consisting of all functions antisymmetric with respect 

to permutations of the n particles. 

As in the bosonic case, we let H = { 1 , 0 , 0 ,... } denote the vacuum and define the number of 
particles operator by (AAT)^”) = for all 'k = Moreover, for all / G L^(M^), we 

introduce the creation and annihilation operators a*{f) and a(/) by setting 

1 ” 

(a*(/)^')(’^)(a:i,... ,x„) = ^ ... ,Xj_i,Xj+i,... ,x„) , 

1=1 (6.19) 

(a(/)T)("-)(xi ,... ,Xn) = y/n + l j dx f{x)y^'^~^^\x,xi,... ,Xn) , 

for T = £ -F. It is easy to see that a*{f) is the adjoint of a{f) and that fermionic creation 

and annihilation operators satisfy the canonical anticommutation relations 

{a*if), a*ig)} = {a(/), 0 ( 5 )} = 0 , {a*(/), 0 ( 5 )} = {g, f) 

for all f,g£ L^(M^). (The anticommutator is defined as {A,B} = AB + BA for any two operators 
A and B.) An important consequence of the canonical anticommutation relations is that a(/) and 
a*if) are bounded operators (unlike in the bosonic case), since 

\\aif)n^ = (T,a*(/)a(/)T) = -(T, a(/)a*(/)T) + ||/||i < ||/||i VT G .A . (6.20) 

It is simple to check that ||a*(/)|| = ||a(/)|| = ll/lb- 

As in the bosonic case, it is useful to introduce operator-valued distributions a* and Ox, which 
allow us to write 

aif) = J dx fix)ax , a*if) = J dx fix)ax . ( 6 . 21 ) 

The second quantization of operators on the Fock space is defined exactly as for the bosonic case; 
see Section [3] 

Fermionic Bogoliubov transformations. It turns out that A^-particle Slater determinants can be 
obtained in the Fock space F by applying appropriate Bogoliubov transformations on the vacuum 
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vector 0. To explain this point in more detail, let us define fermionic Bogoliubov transformations, 
similarly as we did in Section 0] for the bosonic case. 

For f,gG we define 

= a{f) + a*{9)- 

Observe that 

A*{f,g)=A{Jif,g)) (6.22) 

where the operator J is defined by J[f-,g) = ( 5 ,/) (as in (jl3])). In terms of these operators, the 
canonical anticommutation relations assume the form 

{^(/i, 5 i),^*(/ 2 , 5 ' 2 )} = {a(/i),a*(/2)} + {a*(51), 0(52)} 

= (/i,/ 2 ) + ( 52 , 51 ) = (/i,/ 2 ) + ( 51 , 52 ) (6.23) 

= ((/l,5l), (/2,52))L2eL2- 

A fermionic Bogoliubov transformation is a linear map v : L^(M^) 0 L^(M^) —)• L^(M^) © L^(M^) 
preserving (j6.22p and (|6.23p in the sense that these properties continue to hold for the new field 
operators B{f,g) := A{i'{f, g)). It is easy to see that is a fermionic Bogoliubov transformation 
if and only if v is unitary and = vj. Equivalently, a linear map v : L^(M^) © L^(R^) —)• 
L2(M3)©L2(M3) is a fermionic Bogoliubov transformation if and only if it has the form 


u V \ 

V u J 

where u,v : L^(M3) linear maps with 

* I * 1 ^— I *— 

U U + V V = 1 , U V + V u = U. 



(6.24) 


(6.25) 


As in the bosonic case, a fermionic Bogoliubov transformation ly is said to be implementable on B if 
there exists a unitary operator Ry : F ^ F such that 


RlA{f,g)Ry = A{u{f,g)) (6.26) 

for all f,g G L^(]R3). By the Shale-Stinespring condition, the fermionic Bogoliubov transformation 
ly is implementable if and only if u is a Hilbert-Schmidt operator (see e. g. [ 6 OI Theorem 9.5] or [57]]). 

It is useful to see how the unitary Ry transforms the fermionic operators and the corresponding 
operator-valued distributions. We have: 

R*ya{f )Ry = R*yA{f, 0)Ry = a{uf) + a*{vj) (6.27) 

for all / G L^(M3). Setting a{f) = f dxaxf(x), Eq. ()6.27p is equivalent to 

Rla^Ry = a{ux) + a*{vx) , (6.28) 

Here we introduced the notations Ux{y) ■= u{y;x) and Vx{y) = v{y;x), where u{x;y) and v{x;y) are 
the integral kernels of the operators u and v. 

Quasi-free states. For any implementable Bogoliubov transformation v : L^(M3) © L^(m3) ^ 
L^(M3) 0 L^(M3), the Fock space vector Ryid describes a quasi-free state; all correlation functions 
of Ryit can be computed through one-particle correlations by means of Wick’s rule. In this sense a 
quasi-free state T G T" is completely determined by its one-particle reduced density with the 
kernel 
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and by its pairing density , defined by 

a%\x,y) = {^,ayax'^) 

Let wtv denote an orthogonal projection on with Trcj^v = ^■ Then we find an orthonormal 

system {/jjjLi in such that ujn = We define = 1 — ujn and vjq = 

^^1 \ fj){fj\- Observe that u*j^un + v’^vn = 1 and uyvjq = v’^un = 0. Hence (|6.25|) is satisfied 
and 

f UN VN \ 

VN = \ - 

\VN UN J 

defines an implementable (since Ttu^utv = Ttwtv = N < oo) Bogoliubov transformation. This 
implies that Ruj^^ is a quasi-free state. Its one-particle reduced density has the integral kernel 

'r^N{x;y) = {Kj^^,a*yaxRuj^n) 

— (H, [o {uN^y') T (^{vN,y)][(^{uN^x') T Q. (ujv,x)]fl) 

= (H,{a(hAr,y),a*(i;Ar,j:)}H) = {v*^VN)ix]y) 

= ujN{x]y) 

The pairing density of R^^^ is given by 

= (H, [a{uN,y) + a*{vN,y)][a{uN,x) + a*{vN,x)]^) 

— (H, {(i(?r/vy), Q, 

= (vNUN)ix-,y) = 0 


Having the same one-particle reduced density and pairing density as the A^-particle Slater de¬ 
terminant ^/’siater(3^i) ■ ■ ■ ,xn) = det(/j(a:j)), we conclude that up to a trivial phase Rpj^^ = 

{0,0,... ,'i/’Slater,0,0. . . }. 

Alternatively the unitary operator R^^ can be introduced as a particle-hole transformation on 
R. As above, let con = orthonormal system Moreover, let us complete 

{fj}j^n to an orthonormal basis of L^(]R^). Then, Eq. (I6.27p implies that 


Rur^aim* = Rl^aim,* = Rla{fi)R 


'‘UN 


a{fi) ifi> N 

a*{fi) ifi < N 


(6.29) 


Being unitary, Ri,^^ defines new fermionic creation and annihilation operators b*{f) = Run^^* imiN 
and b{f) = Ri,^a{f)R*^, whose vacuum is the Slater determinant Run^- 

Dynamics of quasi-free states. As in the bosonic case, we define the Fock space Hamilton operator 

Rn = J dxVxalVxttx J dxdyV{x - y)a*ayayax ■ (6.30) 

Notice that on the A^-particle sector oIJ-,TLn coincides with the Hamilton operator appearing on the 
r. h. s. of the Schrodinger equation (I6.12p . The next theorem describes the time-evolution generated 
by (jb.OOD on approximate Slater determinants. 
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(6.31) 


Theorem 6.1. Let V G L^(R^) such that 

jdp\V{p)\{l + IpP) < oo. 

Let ujn be a sequence of orthogonal projections on L^(R^) with Truj^ = N and 

Tr\[x,ijJN]\ < KNe, Tr|[eV, cjat]! < i^iVe , (6.32) 

for some constant K > 0. Let be the fermionic Bogoliubov transformation associated with the 
projection un (as constructed in the last paragraph). 

Let fN be a sequence in F with < C and denote by the one-particle reduced density 

associated with the Fock space vector 

Then there exist constants C,c > 0 such that 

^ C'exp(cexpc|t|) (6.33) 

for all i G R. Here uJM,t is the solution of the time-dependent Hartree-Fock equation 

iedtOJN,t = [-e^A -\- V * pt - Xt, ujt] (6.34) 

with initial datauiNfi = ojn (recall that pt{x) = N~^ujn{x]x) and Xt{x\y) = N~^V{x — y)oJM{x\y)). 
Under the additional assumptions that {f,N, ■hf'^^. n) < C and dT{ujN)f,N = 0, we also have the trace- 
norm bound 

~ ^N,t\ < exp(cexpc|t|) . (6.35) 

Remarks. 

• The r. h. s. of (|6.33l) and (|6.35l) should be comoared with the size of IItSJiIhs ~ and of 

Tr 7 ^j = A. In this sense, (I6.33p and (|6.35p prove that the Hartree-Fock equation gives a good 
approximation to the many-body evolution in the fermionic mean-field regime. 

• The vector describes initial deviations from the Slater determinant. Taking = Ll, the 
initial data is exactly the A-particle Slater determinant with reduced density lon- As long as 
(^ 7 V,AA^ 7 v) < C uniformly in A, the initial data Rufff,N is still close to a Slater determinant, at 
least in the sense of reduced densities (choosing t = 0 in (|6.33|1 we observe that ||7^g —watUhs < 
C). 

• It is simple to extend (I6.33p to & F with {f,N,-hff,N) < A“ for any a < 1 (if a = 1, 
the deviations from the Slater determinant are as big as the Slater determinant). Under this 
assumption, ()6.33D is replaced by the weaker bound 

- WAT,dins < CA"/^ exp(cexp(c|t|)). (6.36) 

• It is possible to extend (I6.36P to arbitrary A-particle initial data ipN with reduced density 
close to the one of a Slater determinant in the trace norm topology. Let a ;at be a sequence of 
orthogonal projections on L^(R^) with Trcu^r = A, satisfying (|6.32p . Let ipN be a sequence of 
normalized functions in L^(R^'^) with one-particle reduced density 7 ^^ satisfying 

TrlyJJ^-wjvl < CA" (6.37) 
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for some a < 1. Set = R*^{0, 0,... , i/jn, 0,... } and observe that, using (|6.28l) . 

(Cat, AA^at) = {-ipN, (AA - 2dr(a;Ar) + N)'ijjN) . (6.38) 

Here, we identified V’iv with the Fock space vector {0,0,..., '4’n,0, ■ ■ ■} and we used the identity 

Ru^MRI^ = RI^MR^^ =M- 2dT{ujM) + N (6.39) 

which follows from (16.291) . applied io M = a*(/j)“(/?')• From (|6.38p we find 

iiNi-N'iN) = 2 (V’Ar,dr(l -LOn)'iPn) 

= 2 Tr 7 ^^(l - ujn) = 2Tr ( 7 ^^ - a;Ar)(l - ujn) < 2Tr | 7 ^^ - ujn\ 

where we used that ujn is an orthogonal projection. Hence ()6.37p implies {Cn,-^Cn) < CN°‘ 
for some a < 1, and (I6.36p holds. 

• Theorem 16.II can be extended to prove the convergence of the higher order correlation functions 
7 ^^^ to dehned by the integral kernel 

k 

io^^l{xi,...,Xk;yi,...,yk) = (^nYluiNA^j'^ynU)) ■ ( 6 - 40 ) 

■n-eSk i = l 

(This is just the /c-particle reduced density as it looks like for a quasifree state due to the Wick 
theorem.) More details about this point can be found in [151 Theorem 2.2]. 

• It is well-known that the exchange term is subleading with respect to the direct one. In 
particular, for the class of potentials we consider in Theorem 16.11 dropping the exchange term 
does not change the estimates (I6.33p . In fact, recalling that Xt has the integral kernel Xt{x] y) = 
N~^V{x — y)(jJN,t{x]y), we write 

X* = 1 / dpV{p)e^P-^u;N,te-^^'^ (6.41) 

and we find 

Tt\[Xt,0JN,t]\ <^J dp\V{p)\Tv\[e^P--UN,te-^^-^,u;N,t]\ 

<§/ dp\Vip)\Tr\[e^P-^,ujN,t]\ (6-42) 

<^/ dp\V{p)\\p\TT:\[x,ujN,t]\ 

where in the last step we used that Tr |wat,*]] < |p|Tr |[x,a;Ar,i]|. As we shall see later, 
the bounds (|6.32l) can be propagated along the flow of the Hartree-Fock equation; the result is 
Tr I[x,a;Ar,i]| < KNeexp{c\t\), and hence 

TT\[Xt,ujN,t]\ <Ceexp{c\t\) j dp\V{p)\\p\ . (6.43) 

Roughly speaking, (16.431) implies that the contribution of the exchange term to the solution of 
the Hartree-Fock equation (I6.34p is of order e = (over times of order one), and hence 

comparable with (or of smaller order than) the r. h.s. of (|6.33p and (|6.35p . The situation is 
expected to change in the presence of unbounded potentials, e. g. for a Coulomb potential; in 
this case, the analogy with the study of the ground state energy of large atoms [ElET] suggests 
that neglecting the exchange term will deteriorate the error estimates. 
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Sketch of the proof of Theorem \ 6.1[ The proof follows a strategy conceptually similar to the 
coherent state technique discussed in the bosonic case in Section [31 For simplicity, we shall only 
discuss the proof of convergence in Hilbert-Schmidt norm, that is of the inequality (I6.33p . Let 
'^N,t = be the evolution in Fock space of the initial state The kernel of the 

one-particle reduced density 7 ^ ^ is 

y) = R.^^n) • (6.44) 

The idea is to compare the quantum evolution with the Hartree-Fock evolution, de¬ 

scribed in Fock space by Ri,ff^f,N- This is done by introducing the fluctuation dynamics UM{t', s) := 

Using the shorthand notation ut^xiu) = UN,tiy,x) and vt^xiu) = VN,tiy,x), 
we can rewrite the kernel of as 

, Rl^ _^a*R„j^_^Rl^ _^axRuNtUN{t-, 0 )f,N) 

= {UN{t-,0)^N, [a*{ut,y) + a{vt,y)][a{ut,x) + a*{vt,x)]h{N{t-,0)f,N) 

= u}N,tix;y) (6.45) 

-h {UNit] 0 )(,N, [a*iut,y)a{ut,x) + a*{ut,y)a*{vt,x) 

+ a{vt,y)a{ut,x) - a*{vt,x)a{vt,y)]UN{fl fl)fN) ■ 


Here we used the transformation property (|6.28p and the fact that {a*{vt^x),<i{vt^y)} = ujN,t{x'-,y)- 
For any Hilbert-Schmidt operator O on we find 


TrO( 7 /^,t “ ^N,t) = [dT{uN,tOuN,t) - dr{vN,tOvN,t)] ^Af(i; 0 )^Ar) 




dxdy 0(x; y){UN{t] 0)fN, a{vt,y)a{ut,x)hlN{t; 0)f,N) + c.c. 


(6.46) 


Consider the first two terms in (I6.46p . Since for any bounded self-adjoint operator A on L^(M^) 
we have the bound ±dF(H) < ||H||AA, where ||H|| is the operator norm of A, we find, using that 
||wAr,t|| < 1 and ||uAr,t|| < 1 , 


{UN{t-,0)f,isf, [dr{uN,tOuN,t) - dT{vN,tOvN,t)]h(Nit; 0 )f,N) < ‘^\\O\\{UN{t-,0)f,isf,N'UN{t-,0)f,isf) . 

(6.47) 

Consider now the second line of (I6.46p . In order to estimate it, we shall use the following bound. 
Let H be a Hilbert-Schmidt operator on L^(M^), with kernel A{x',y). Then 


J dxdy A{x-,y)axay{M + 1) 


■ 1/2 




< 


dy 


a{A{-;y))ay{M +1) 


< j dy\\A{-;y)\\ 2 \\ay{J\f + 1) U2^|| 

< PIIhs j rfy||ay(AA-hl)”U2^|| 

= ||H||hs|KU2(a7 + 1)-U2^|| < IIHIIhsIIV'II ■ 


(6.48) 


where in the second line we used the boundedness of the fermionic creation/annihilation operators, 
Eq. (I6.20p . We find 


(AA-l- 1 ) J dxdy A{x-,y)a*a*fj = J dxdy A{x;y)aya*{N' + 3) 


< 


HSi 


. (6.49) 
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The estimates (|6.48l) . (16.491) can be used to bound the second line of (16.4611 . We get 


I dxdy 0(x; y){UN{t, 0)^N,a{vt,y)a{ut,x)ldN{t; 0 )^n) + c.c. 

< C\\uN,tOvN,t\\KS ||(AA+ l)^/‘^UN{t]0)^N\\ ■ 
Since ||ttiv,tOuAr,i||HS < ||0||hs and ||0|| < ||0||hs, (|6.47|) and (I6.50p imply 

Ti’0(7^t — a;Ar,i) — C'||0||Hs(^Ar(i; 0).^Ar, (W + 1)7/Ar(t; 0),^7v) • 


(6.50) 


(6.51) 


In particular, choosing O = 7 )^') — we obtain 

llo'w,! “< C{UN{t;0)^N,iJ^ + i)^N{t;0)^N) ■ (6.52) 

Hence, to conclude the proof, we need to control the growth of the expectation of the number of 
particles operator with respect to the fluctuation dynamics 0). This is the content of the next 
proposition. 


Proposition 6.2. Under the same assumptions as in Theorem 16.11 there exist constants C,c > 0 
such that 


(6/Ar(t;0)^Ar, (AA+l)6/Ar(t;0)^Ar) < Cexp(cexp c|t|)(^ at, (AA + l)^Ar). (6.53) 


Controlling the growth of the fluctuations. Here we present the main steps in the proof of the key 
bound Eq. (j6.53ll . It is based on a Gronwall-type argument. Using the identity (I6.39p we write 


= [Af + N- 2dr(a;Ar,i)]e-'^^'/^i?,^^Ar). 


(6.54) 


The expectation of (Af+N) is independent of time; these terms disappear when we take the derivative. 
Hence 


0)f,N,AfUNit; 0)f,N) 

= {dT{i£dtOJN,t) - [RN,dT{ojN,t)])e-^'^^^/^R,^fN) 

= ,RvNy{dT{i£dtUJN,t) — \R-N ,dT{wN,ty\) RuNyhlN{t',^)f.N) 

= -‘^^N{t]^)fN,RvNy{dT{pt*V - Xt) - [VAr,(ir(wAr,t)])i?i/jv_tZ7Ar(t;0)^Ar) 


where 



dxdy V{x — y)a*a*ayax 


is the many-body interaction. In the last step we used that the contribution of the kinetic en¬ 
ergy cancels, as follows from the identity dT{[A,B]) = [dr(H), dr(i?)]. Applying the Bogoliubov 
transformation R^y^.^ ^ and reorganizing all terms in normal order we end up with 


ie— (U]si{fl U)f^N,AfUN{t] 0)^Ar^ 

Ai f 

= — Im / dxdy V{x — y) 

X ^(uN{t;0)f,N,a*{ut,x)a{vt,y)a{ut,y)a{ut,x)h(Nit-,0)f,N^ 

-h (uNifl 0)f,N, a{vt,x)a‘{vt,y)a{ut,y)a{upx)UN{t] 

+ (uN{t]t))^N,a*{ut,y)a*{vt,y)a*{vpx)a{vt,x)UN{flQ)iN^Y 
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Consider for example the second term on the r. h. s. of (|6.56l) . Expanding the potential in Fonrier 
space, we write 

^ j dxdyV{x - y)(uNit]0)^N,a{vt,x)a(vt,y)a{ut,y)a{ut,x)l^Nit‘,0)^N) 

= ^ J dpV{p) 

X j dridr2dr3dr4{vte~^P'^ut)iri;r2){vte'-P'^ut){n-, r4)(uN{t] 0)^Ar, ^r2 0)^A^ ^ • 

Inserting 1 = (A^ + 5)^/^(Ar + 5)”^/^ and applying the Cauchy-Schwarz inequality, we get 

J dxdyV{x - y)(uN{t;0)^N,a(yt,x)a{vt,y)a{ut,y)a{ut,x)idN{t;0)^N) 

^ j dp\V{p)\\\iAf + 5)^/‘^UNit-,0)^N\\ 

X (AA + 5)"^/^y dridr2ivN,te~"^'^UN,t){ri]r2)ariar2 

X J dr3dr4{vN,te'-P'^UN,t)ir3;ri)ar3ariUNit;0)^N 
= jjj dp\V{p)\\\{M+ 5)^/‘^UN{t;0)CN\\ 

X j dridr2ivN,t[(^~'‘^'''i'^N,t]){ri;r2)ariar2 

X j dr3dr4{vN,t[^"^'''^UN,t]){r3;r4)ar:iar^{J\f + l)~^^‘^UN{t-,0)^N 

where in the last step we used that Ma{f) = a{f){M — 1) and the orthogonality vtUt = 0. Applying 
the estimate ()6.48p twice and using that ||nAr^i|| < 1 we get: 


1 

iV 

< 


y J dxdyV{x - y)(uNit;0)^N,a{vt,x)a(vt,y)a{ut,y)a{ut,x)idNit]0)^N) 

< ^ J dp\V{p)\\\[e''P''^ ,u:N,t]\\Rs(^^N{t',0)^N,J^idNit-,0)^N 


(6.57) 


Using the naive bound ||[e*^'*, WAr,f]||Hs — 2|| WAr,t] ||tr < dA" would not be enough here because of 
the factor s on the 1. h. s. of (|6.56li . Instead, we have to propagate the semiclassical estimates (I6.32p 
along the solution of the Hartree-Fock equation. This is the content of Proposition 16.31 below, which 
tells us that ||[e®^'^,a;Ar^ 4 ]||tr < CNe{l + |p|) exp(c|t|) for all times t G R. Inserting this estimate in 
(I6.57p . we conclude that 

y j dxdyV{x - y)(uNit;0)^N,a{vt,x)a(yt,y)a{ut,y)a{ut,x)idNit]0)^N) 

< C£ex.p{c\t\){L{N{t; 0)^N,dd'ldN{t-, 0 )^n)- 


The other terms in (I6.56P can be bounded in a similar way. The only difference is that instead 
of (I6.48P one has to use the estimate 


dr I dr 2 A (ri; r 2 ) aj, ar^ y:> 


< 


tr 


Vv? G A 


(6.58) 
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valid for any trace-class operator A on L^(M^) with kernel A(ri;r 2 ). The details can be found in 
m- We conclude that 


A 

dt 


UN{t-,o)n,{Af + i)UNit]0)n) <Cexp{c\t\)(uNit-,o)Q,{M + i)UN{t-,o)n). 


(6.59) 


Gronwall’s lemma implies the desired estimate (j6.53p . 

Propagation of the semiclassical structure. An important ingredient that we used above to control 
the growth of the expectation of the number of particles with respect to the fluctuation dynamics 
is the propagation of the semiclassical commutator bounds along the solution of the Hartree-Fock 
equation. 


Proposition 6.3. LetV satisfy li6.Sl\) . Letups he a sequence of orthogonal projections with Truipf = 
N, such that 

TrI[x, Wat]! < ATiVe and rr|[eV,a;Ar]| < (6.60) 


for some constant K > 0. Let ujN,t be the solution of the Hartree-Fock equation \6.34\) with initial 
data ujpffi = ojn- Then there exist constants C,c > 0 such that 


Tr|[x,a;Ar^t]| < CiVeexp(c|t|) and rr|[eV,a;Ar^t]| < CiVeexp(c|t|) 


for all t G M. 

Sketch of the proof of Proposition \ 6.A Let hHpit) = —e^A V * pt — Xt he the Hartree-Fock 
Hamiltonian. We compute 


ie-^[x,ujN,t\ 


[x,[hHF{t),U]N,t]\ = [^N,t,[hHF{t),x]] + [hHF{t),[x,U]N,t]]- 


The last term can be eliminated by conjugating [xjWtv,*] with the two-parameter group of unitary 
transformations W (t; s) defined by 


ie-^W (t; s) = hHF{t)W{t; s) with W (s; s) = 1 for all s G M. 

We have 

iejW*{t-,0)[x,iON,t]W{t;0) = W*{t-,0)[iON,t,[hHF{t),x]]W{t;0) 

= W*{t-,0)([ujN,t,-2eW] - [ujN,t, [W,x]])w(t;0) 

where we used the identities [—e^A, x] = —2e^V and [pt *V,x\ = 0. Therefore we get the Duhamel- 
type formula 

W*{t,0)[x,u;N,t]Wit-,0) = [x,a;7v,o] - ^ dsW*{s-,0){[uN,s,2e^X] + [A„ x]]) W(s; 0). 

This implies that 

||[x,a;Ar,i]||tr < || [x, ||tr + ^ y ds (^||[wAr,s,2e^V]||tr + ||[wAr,s, [X5,x]]||try (6.61) 
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To control the second term we use (j6.41l) . Since ||wAr,s|| < 1 we find 


||[u;^,„[X„x]]||tr <^l dq\Viq)\ ||[cn^,„ 

<^/ dq\V{q)\\\[ujN,s,e^^'^[ujN,s,x]e~^P'^]\\tr 
j dq\V{q)\ ||[a;Ar,^,x]||tr 

c 

— ^ II 3^] ||tr- 

Inserting (16.6211 in (|6.6ip we get 

II [x, ||tr ^ II [x, CJjVjO] ||tr T C ds (11 [u;n,s, e V] 11 tr + 11 [uJN,s , x] 11 tr 


To control [a;Ar,s,eV] we start by writing 


(6.62) 


(6.63) 


is—[sV,ujN,t] = [£V,[hHF{t),u:N,t]] 

= [^N,t, [hHF{t),£'^]] + [hHF{t), [eVjWAT^t]] 

= [hHF{t), [£V,ujN,t]] + [u}N,t, [pt * l^,eV]] - [ujN,t, [A't,eV]]. 

As before, the first term can be eliminated by conjugation with the unitary operator W{t;0). We 
find 


||[eV,a;Ar,t]||tr < ,^N,o]\\tr + ^ j ds ( ||[wAr,s, [ps * y,eV]]||tr + ||[wAr,s, [^s,eV]]||tr )• (6.64) 


The first term in the integral can be controlled by 


||[wAr,s,[y *Ps,eV]]||tr = ellfwAT,*, vy */ 9 s]||tr 

<e [ dg|y(g)||g| |p,((?)| || ||tr 

(6.65) 

<£j dq\V{q)\\qf\ps{q)\\\[u;N,s,x]\\tr 
— f^^ll lltr 

where we used the bound ||/5s||oo < lids 111 = 1 and the assumption on the interaction potential. 
Finally, consider the second term in the integral on the r. h. s. of (j6.64p . Again using (I6.4ip . we write 

||[a;^,„[A„£V]]||tr <^l dq\V{q)\ ||[a;^,„ [y‘^-a;jv,.e-*''-,eV]]||tr 

< I / dq\V{q)\ ||[y'?-a;^,,e-*''-,£V]||tr (6.66) 

j ^^9l^(9)l Ill^iV.sy'^llltr 

where we used the identity 
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Inserting the estimates (|6.65l) and (I 6 . 66 P into (|6.64l) . we get 

||[eV,u;Ar,t]||tr < ||[eV,WAr^oJUtr + C J ds (^||[a;Ar,^,x]||tr+ iV"^'^^||[wA^,^,eV]||tr)• 

Combining this ineqnality with (I6.63|) , using the assumptions on the initial data and applying Gron- 
wall’s lemma, we obtain (jb.GOp . 

7 Dynamics of quasi-free mixed states 

In Section [ 6 ] we discussed the evolution of initial data approximating Slater determinants. Slater 
determinants are relevant at zero temperature because they provide (or at least they are expected to 
provide) a good approximation to the fermionic ground state of Hamilton operators like () 6 .ip in the 
mean-field limit. At positive temperature, equilibrium states are mixed; in the mean-field regime, 
they are expected to be approximately quasi-free mixed states. 

Mixed states. We introduce the short-hand notation i) = L^(]R^). We shall denote by the 
fermionic Fock space built over t), that is A general fermionic state is represented 

by a density matrix on J^(l)). A density matrix is a non-negative trace class operator p:J^(f))— 
with Trp = I. Notice that the state described by the density matrix p is pure if and only if p is 
a rank-one orthogonal projection onto a ip ^ he. p = \ip){ip\. Otherwise the state is called a 

mixed state. In general, 

/5 = ^An|V’n)(V’n| (7.1) 

n 

where \n > 0, {ipn} is an orthonormal family in 7^(1}), and = 1- Physically, p describes an 

incoherent mixture of pure states and is the probability that the system is in the state ipn- The 
expectation of an arbitrary operator A on in the mixed state with density matrix p is given by 

TrAp = Xnjlpn, Alpn). 

n 

Given the density matrix dm) we define the operator k : ^ by 

n 

where £„ G C is a sequence satisfying |e„p = A^ and {(pn) is an orthonormal family in J^(l)). Glearly, 

'kR* = p. 

Of course such a decomposition of p is far from being unique and later we shall choose a convenient 
one. Since p is trace class it follows that R G the set of Hilbert-Schmidt operators on 7^(1)). 

Next, we observe that £^(J^(f})) ~ 7^(f)) <8* F{^). This isomorphism is induced by the map 
|V')(())| Ip <Si (p, extended by linearity to the whole space i2^(J^(f))). The mixed state with density 
matrix dZU) can be described on 7^(1)) 0 F{\]) by the vector 

K = ® (Pn- 

n 

The expectation of the operator A on F{\)) in the state k G 7^(1)) ® -T(f}) is given by 

Tr Ap = Tr ARR* = {k, {A (g) I)K)jr([,)^jr((,) . (7.2) 
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The doubled Fock space (8) is isomorphic to T'(() © fj) (see [23] or any book on mathe¬ 
matical quantum field theory). The unitary map U that implements this isomorphism is known as 
the exponential law and is defined by the relations 


and 


U[a{f)®l]U* = a{f(S0)=:ai{f) 

U [(-1)^ © a(/)] U* = o(0 © /) =: a,(/) 

for all / G f), where a^if), a = I, r, are called the left and right representations of a(/), respectively. 
By hermitian conjugation, we also find 


[/[a*(/)©!] U* = a*{f®0) =:<(/) 
U [(-1)^ © a*{f)] U* = a*(0 © /) =: <(/) 


(7.4) 


where a* (/), a = I, r, are the left and right representations of a*{f). Notice that the presence of the 
operator (—1)'^ on the second line of ()7.3p and (17.41) guarantees that creation operators on the space 
T'(f)© fi) satisfy the canonical anticommutation relations (and in particular that aj{f) anticommutes 
with at{g), for all /, <7 G fl). 

It is convenient to introduce the left and right representations of the operator-valued distributions 
ax, a* by the relations 


aiif) = ^ 

j dx ax^ifix), 

arif) = ^ 

j dxOx^rfix). 


[ dx allfix), 

a*rif) = ^ 

[ dx al J{x), 


(7.5) 


for all / G f). We also define the left and right representations of the second quantization of operators 
on f) by 


U [dr(0) ©!][/* = dr{0 © 0) =: dViiO), 

[/ [1 © dr(0)] U* = dr(0 © O) =-. dVriO) . 

The left and right representations of dr(0) can be written in terms of the operator-valued distribu¬ 
tions as 

dTi{0) = j dxdy 0{x‘,y)a* iay^i, dTr{0) = J dxdy 0{x]y)a*^^ay^r- (7.6) 

According to (|7.2p . the expectation of the observable A in the mixed state described by the vector 
-0 = Uk G © pi) is given by 

{k, {A © l)K)_7r((,)0_7r(f,) = {Uk,U{A® 1) [/* I/k) = {if, U {A ® 1) [/* V') . 

In particular, the expectation of the second quantization (ir(0) of a one-particle operator O on f) is 
given by 

{ii,U{dT{0) © = {'if,dTi{0)'if) = j dxdy 0{x-,y){if,al^iayj'if) = TrOy^^^ 
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where we defined the one-particle reduced density 7 ^^^ associated with ip G as the non¬ 

negative trace class operator on f) having the integral kernel 

7^^^ {x;y) = {'ip, a* ia^^iip). (7.7) 

We also define the pairing density associated with ip ^ as the Hilbert-Schmidt operator 

on f} with the kernel 

a.^{x\y) = {ip,ay^iax,iip). (7.8) 

The above construction allows us to represent mixed states as vectors in the “larger” Fock space 
T'(f)©f)). This idea is well-known in quantum statistical mechanics and takes the name of purification. 

Time evolution of mixed states. The time evolution of the density matrix p is given by 

Pi = 

where 7 ^tv is the second quantized Hamiltonian given by (|6.30p . Accordingly, we define the time 
evolution of k € T'(f)) ® T'(f)) by Kt = ® g*'K]v7ej Let ipt = U denote the vector in 

T'(f) © {)) describing the mixed state with density matrix pt. Then 

= UKt = 


where the Liouvillian Ttv is defined by 

£jv = C/(7 ^tv® 1 - 1©7^tv)C7*- 
A more explicit expression follows from (j7.3p and (17.5^ : 

Cn = J dxV^al iVxa^^i + ^ j dxdyVix - y)al ia*y iay^iaa,^i 


— e 


dx^xO>x,r^^^oc,r 2jy ^ 


(7.9) 


Hence, the expectation of an arbitrary operator A on T'(f)) in the evolved mixed state is given by 
TrApt = {ipt,U{A(^l)U*ipt) = {ip,e^^^*/^U{A®l)U*e-^^^*/^ip). 

Araki- Wyss representation. In Section [6l where we proved that the time-evolution of Slater 
determinants can be approximated by the Hartree-Fock equation, a crucial ingredient of our analysis 
was the observation that Slater determinants can be written in the Fock space T'(f)) in the form 
where i'n is an appropriate Bogoliubov transformation. In this section, we are interested in 
the evolution of quasi-free mixed states, which, after purification, can be described by vectors in the 
Hilbert space T'(f) © f)). 

Let u : (t)©f})©(f)©f)) —>• (f)©l|)©(f)©li) be an implementable Bogoliubov transformation on the 
doubled one-particle space. Then there exists a unitary map Ri, : J-'(f)©f)) —>■ T'(t)©t)) implementing 
n. It is easy to see that vectors of the form G 7^(1} © f}) are quasi free in the sense that 

higher order correlations can be computed through Wick’s rule. 

In the following we will be interested in quasi-free states in J-'(t) © tj) with average number of 
particles equal to N and with vanishing pairing density. Let a ;at be a non-negative trace class operator 
on f) with 0 < ujn < 1 and Trwjv = A^- Then we define i'n ■ (f) © ()) © (f) © f)) —>■ (f) © f)) © (f) © ()) by 


UN = 


Un Vn 
Vn Un 


(7.10) 
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where 

U„ = = ( “ 7) , {7.11) 

VO UN J \-VN 0 J 

and UN = VI “ ^N, vn = V^Vv- The operators Un, Vn : f) © f) ^ f) © f) satisfy the relations (16.2511 : 
hence i^n defines a Bogoliubov transformation. Since TtV^Vn = 2Tru;7v = 2iV < oo, the Bogoliubov 
transformation un is implement able. Hence ^ .T(f) © f)) describes a quasi-free mixed states 

(here and in the following, we denote simply by H the vacuum in J^(f)© f))). We claim now that 
is exactly the quasi-free mixed state with reduced density lon and with vanishing pairing density. In 
fact, for any / G f), we find 


Ki,(^i{f)RyN = ai{uNf) - a*{vNf), 
KNar{f)RuN = ar{uNf) + a^VN?)- 

Equivalently, 

Rl/j^^x,iRu{^ — Oii{un,x) N,x) ■} 

Ryj^Otx^rRuf^ — (^riuN^x) © O; {un^x')} 


(7.12) 


(7.13) 


with the usual notation UN,x{y) = UN{y',x), VN,x{y) = VN{y',x). Using (17.1311 and the unitarity of 
Ruiq we find, from (17.7D . 


y) - {Ruff^T 0,y^l0‘X,lRuj^^) 

= {n,R:^aliR,,R:^ax,iRu,n) 

= (H, {a*{uy) - ar{vy)){ai{ux) - a*(yx))i^) (7-14) 

— (H, (n^)(z^ (n^,)!!) 

= {v*nvn){x; y) = tONix] y) 


where we used the canonical anticommutation relations and the fact that ai and annihilate the 
vacuum. Also, from (j7.8p . 


(x; y) — {R^j^Vl,ay^i(ix,lRi'j^^) 

= {Q, {ai{uy) - a*{vy){ai{ux) - a*(vx))fl} (7.15) 

= (fl, ai(uy)a*(vx)fl} = 0 

where we used that {ai(uy), a*(v^)} = 0. 

This representation of quasi-free mixed states is a special example of a well known construction 
in quantum statistical mechanics, known as the Araki-Wyss representation [6j. 

Dynamies of quasi-free mixed states. We consider the time evolution of initial quasi-free mixed 
states satisfying certain semiclassical estimates (motivated by the idea that, physically, we are inter¬ 
ested in the evolution of equilibrium states at positive temperature). For such initial data, we show 
that the evolution remains close to a mixed quasi-free state with reduced density evolving according 
to the Hartree-Fock equation (16.611 . The next theorem is taken from m- 

Theorem 7.1. Let V G L^(M^) and assume that 

j dp{l + \p\‘^)\V{p)\ < oo. (7.16) 
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Let ujn be a sequence of operators on \] = L^(M^) with 0 < lon < 1, Truj^ = N and such that 
Tr{l — A)ujn < oo and 


IIK,x]|| HS < CN^/^s, \\[vn,£'^]\\hs<CN^/^s, 

||[niv,x]||i^5 < ||[n^,£V]||^^5 < CN^/^, 


(7.17) 


with vn = y/dJN, UN = y/T^-ojN, for a suitable constant C > 0. Let vn denote the Bogoliubov 
transformations (f7.70D , such that is the quasi-free state on Rif) 0 f)) with one-particle 

reduced density ojn and with vanishing pairing density. 

Let 7 ^^ j be the one-particle reduced density associated with the evolved state 


fjN,t = e 


(7.18) 


where the Liouvillian Cn has been defined in 17.^1 ). 
Hartree-Fock equation 


Let u)N,t be the solution of the time-dependent 


iedtUJN,t = [-e^A -^V * pt- Xt,UN,t] (7.19) 

with the initial data ujn,o = ojn- Then there exist constants C, c> 0 such that 

^5 < C'exp(cexp(c|i|)) and Tr\'y^j^\ — ujN,t\ ^ ex.p{cex.p{c\t\)). (7.20) 

Remarks. 


• Similarly to Theorem 16.11 the convergence towards the Hartree-Fock dynamics can be extended 
to more general initial data than those appearing in (I7.18p . Let Rf = dr;(l) 0 (irr(l) be the 
total number of particles operator in the doubled Fock space and let ^n be a sequence in 
R{i) 0 1 )) with (Cw)< C and such that ^n = x(AA < KN)fN (for a sufficiently large 
constant K > 0 independent of N). In [13] it is shown that the one-particle reduced density of 
Ruff(,N can be approximated by the solution of the Hartree-Fock equation. 

• The semiclassical commutator estimates ()7.17p play the same role as (I6.32p in the analysis of 
the evolution of Slater determinants. There is however an important difference. Eq. (|7.17p 
gives bounds for the Hilbert-Schmidt norm of the commutators, while (I6.32p was expressed 
in terms of their trace norm. In fact, for Slater determinants we do not expect ¥1.17} to be 
correct. In this case the decay of the kernel ojn{x', y) for |x — y| 3> e is quite weak (because ujn 
is a projection and the decay of ww is dictated by the regularity of the Fourier transform in 
the {x — ?/)-direction) and does not suffice, in general, to prove (|7.17p . On the other hand, in 
Theorem 17. 11 ojn does not need to be a projection, and one can expect faster decay of the kernel 
iONix‘,y) in the {x — y)-direction. This is the reason why the assumption ¥1.17} is appropriate 
for the study of the dynamics at positive temperature. 

For instance, a reasonable approximation for the one-particle reduced density of a thermal 
state of the Hamiltonian dSI]) describing trapped interacting fermions is given, in the mean- 
field limit, by the Weyl quantization 

uiN{x\y) = J dvM (7.21) 

where M{x,v) denotes the phase-space density, 

M{x, v) = gT,fj. (v‘^ - cp^F (^)) > (7-22) 
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corresponding to the Fermi-Dirac distribution 


9t,^i{E) ^ ^ , (7.23) 

depending on the temperature T > 0 and on the chemical potential /r > 0. In (|7.22p . ptf is 
the minimizer of the Thomas-Fermi functional (16.711 and the normalization constant c > 0 has 
to be chosen so that 

= 1 - 

For the reduced density wtv, we find 
[wAr,x](x;y) = 

and 

[ujN,£'^]{x-,y) =j dv{V^M) 

Hence, 

||[a;Ar,x]||Hg < e^iV J dxdv\{VyM){x,v)\‘^ < CNs^, 

||[a;Ar,x]||Hs < j dxdv\{VxM){x,v)\'^ < CNe^, 

using the regularity of the Fermi-Dirac distribution ()7.23p (and assuming some regularity of the 
Thomas-Fermi density ptf)- This motivates the assumption (I7.17p for initial data describing 
thermal states at positive temperature (the condition (I7.17P is actually expressed in terms of 
vn = \/djN and un = \/l “ but the same argument is expected to hold also in this case). 

• The proof of Theorem 17.11 follows a strategy conceptually similar to the one of Theorem 16.11 
As for pure states, the rate of the convergence towards the Hartree-Fock evolution can be 
estimated by controlling the growth of the expectation of the number of particles operator 
M = (iF;(l) -|- dFr(l) with respect to the fluctuation dynamics, which is now defined by 

= . (7.24) 

It turns out, however, that new ideas are needed here to control the growth of 

(7/iv(t; 0)D, A/'7/Ar(t;O)D)jF(|,0(,). 

The main difference with respect to Section [6] is that now un = \/l “ and = s/un 
are not orthogonal. To circumvent this problem, one has to use certain cancellations between 
different terms in the generator of the fluctuation dynamics. Moreover, it is important to 
introduce an auxiliary dynamics Z7iv(t; s) which, on the one hand, can be proven to stay close 
to the original fluctuation dynamics, and, on the other hand, only changes the expectation of 
the number of particles in a controllable way. The details of the proof can be found in m- 
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A The Role of Correlations in the Gross-Pitaevskii Energy 


In this appendix we consider bosonic systems in the Gross-Pitaevskii regime, as in Section O We 
show that the energy of Fock space states of the form W{^/N (^)To^jv, for ^ J- with a finite number 
of particles and finite energy describes excitations of the condensate), is to leading order given by 
the Gross-Pitaevskii functional ()5.4|) evaluated on the one-particle wave function (p G This is 

an instructive calculation since it shows that the Bogoliubov transformation Tq converts part of the 
many-body kinetic energy into a contribution to the quartic term in the Gross-Pitaevskii functional. 
Without the Bogoliubov transformation Tq the approximate coherent state W{^/N would have 
a higher energy, given by a functional similar to (1521) but with coupling constant bo = fV(x)dx, 
which is always larger than dvrao- Hence, introducing the Bogoliubov transformation Tq lowers the 
energy by a quantity of order N (while the change of the number of particles is only of order one; 
see ()5.35p b This observation supports the claim that states of the form W(\/N(f)To^]Y provide a 
good approximation for the many-body ground state and that Tq implements the correct correlation 
structure. 

Energy of Bogoliubov states. Consider the Hamilton operator describing a Bose gas in the Gross- 
Pitaevskii regime, trapped by a confining potential I4xt, 


n 


trap 

N 


dx “ 1 “ 2 


dxdy N'^V{N{x — y))a*ayayax. 


Let G T” be such that 

(^fN, (^^f + < C- (A.25) 

uniformly in A^. As in (I5.3ip . we define the Bogoliubov transformation 

To = exp Q j dxdy{kQ{x;y)a*al-ko{x-,y)axay)^ 

with the kernel ko{x', y) = —Nuj{N{x—y))p{x)p{y) . We claim that for sufficiently regular p G 
with ||y5||2 = 1 

(w{^/Np)TQfN,'H^N'’W{y/Np)TQiN) 

= y dx[\Vp{x)\^+ Ve-^t{x)\p{x)\’^+AT:aQ\p{x)\^)+0{'/N) (A.26) 

= A^TgpM + o(Viv). 


Taking p to be the normalized minimizer of TgPj we conclude from [50] that W{y/Np)TQf^j<i has, 
in leading order, the same energy as the ground state of the restriction of to the A^-particle 

sector. 

Sketeh of the proof of l!A.2f)\} . Let us consider the case = H. We compute 

(w{y/Np)TQn,'H%^^W{y/Np)TQdj. (A.27) 

conjugating the Hamiltonian first with the Weyl operator (producing a shift of the creation and 
annihilation operators) and then with the Bogoliubov transformation (which acts on the annihilation 
and creation operators as given by (15.321) 1. At the cost of commutators appearing, all terms can be 
brought into normal order. At the end, all terms with creation and annihilation operators written 
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in normal order vanish when we consider the vacuum expectation. Hence, we find the following 
contributions to (|A.27I) . From the kinetic energy we have 

J dx (^n,T^W{VN^)*V:,alV^a^W{VN^)Ton^ = N J dx\V^{x)\'^ + j dx ||V,,sh,,f . (A.28) 

From the external potential we obtain 


/ 


dxFext(x) (n,T^W{^/N^)*ala,,W{^/Nip)Ton) = N J dx|^(x)|Ve,t(x) + j dxFext(a;)||sha;||2 . 

(A.29) 


Finally, from the interaction, we get 

J dxdyN^V{N{x - y)) (^i},T^W{VN^)*ala*yaya:,W{VN^)Ton 

= \ j dxdyN‘^V{N{x-y))\{ch.y,s\i^)\^ 


+ 


+ 


i J dxdyN^V{N{x - y)) {{sh^,chy)(p{x)ip{y) + c.c.) 
N f 

— J dxdyN^V{N{x-y))\ip{x)\‘^\ip{y)\‘^ . 


(A.30) 


To evaluate the terms on the r.h. s. of (|A.28I1 . (jA.29p and (|A.3nil it is useful (like in Section [5]) 
to think of sh^iy) ^ ko{y;x) and ch^iy) ^ S{x — y); more precisely recall that, in contrast to the 
regular kernels of higher powers of fco, ko{x; y) itself is singular as (|x — y\ + 1/N)~^ for N ^ oo. We 
then find (in the last step using Hardy’s inequality) 

llshxlli^ j dy\ko{y,x)f = \(p{x)\‘^ j dy\ip{y)f \Nu}{N{x - y))\^ 


< 


C\ip{x)\‘^ J 


dy 


|x - ?/P 




This implies that 


dx v;xt(3;)||sha;||2 < C, 


(A.31) 


i. e. the contribution is of order one (for sufficiently regular ip G L^(M^)). Hence, this term can be 
neglected. 

As for the second term on the r. h. s. of ()A.28p . we notice that 

y dx||Vj,sha;|p = j dx{sh^,{-Aj;)sh.^) 

- - j dxdyNu{N{x - y))ip{x)ip{y)A^ [Nu; {N (x - y)) ip{x)(p{y)] 

= J dxdyN^uj{N{x - y))\(p{x)f\ip{y)\^{-Aaj){N{x - y)) + 0{y/N) 

because contributions arising when one or two derivatives act on (p are of smaller order. The key 
observation now is the fact that l — ui satisfies the zero-energy scattering equation (|5.2p , which implies 
that 

{-Auj){N{x-y)) = ^V{N{x-y)){l-u}{N{x-y))) 
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and thus that 


J dx\\Va;sKf = ^ j dxdyN'^V{N{x-y))Lo{N{x-y)){l-uj{N{x-y)))\ip{x)\^\Lp{y)\^+ 0{^/N). 

(A.32) 

Thus, from (|A.28h . (|A.29p and (|A.30p . we conclude (using again the approximations ch.y{z) ~ 
6{z — y) and shj,.( 2 ;) ~ ko{z-^ x)) that 


= TV y dx\S/<f{x)\^ + N J dxVe^tix)\ip{x)\'^ 

N r 

+ ^ j dxdy N^V {N {x - y)) uj {N {x - y)) {1 - ui {N {x - y))) \(p{x)\‘^\(p{y)\‘^ 
N r 

+ Y J dxdyN^V{N{x-y))oj{N{x-y)f\(p{x)\^\(p{y)\‘^ 
j dxdyN^V{N{x-y))u}{N{x-y))\ip{x)\^\(p{y)\^ 

+ y/ dxdyN^V{N{x-y))\ip{x)\^\(p{y)\'^ + 0{y/N)). 

Combining all terms proportional to the interaction potential, we find 


{W{y/Nip)TQQ,'H%'''^ W{y/Nip)TQ9) 


= N 


dx (|V(^(x)p + Ve^tix)\ipix)\'^) 


+ - 


dxdyN^V (N {x - y)) f{N{x - y))|(/?(x)p|(/?(?/)p 


+ 0{^fN). 


Since f dx V(x)f(x) = Svrao, we obtain in the limit N oo (using again some regularity of ip) 


(W(^/Np)Tofl,n^j^^PW(v^p)Ton) = N j dx(|V¥j(x)|2 + Kxt(x) + 47rao|99(x)n 


+ 0{^/N), 


which proves (IA.26p . This computation can be easily extended (estimating all normally ordered 
terms emerging from pA.27p by (jA.25p i to arbitrary satisfying (IA.25h . 
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